MATLAB 图像处理 头颅模型ct图像重建 |
您所在的位置:网站首页 › 头模型图片家里 › MATLAB 图像处理 头颅模型ct图像重建 |
主要任务: 1.建立Shepp-Logan模型,使用灰度叠加 2.生成投影数据 3.反投影重建 4.滤波反投影重建 重建过程会产生星状伪迹,因此需要先修正再投影,不使用Radon函数,原理不讲,直接放代码,模型建立根据原函数源码改编,查看函数源码可在命令行窗口输入: >>open 函数 其余过程不是我自己写的,但是时间过久原网址找不到了 %Shepp-Logan模型 n=256; q=zeros(n); xax=((0:n-1)-(n-1)/2)/((n-1)/2); xg =repmat(xax,n,1); % x坐标,y坐标旋转90度,形成块矩阵 % ρ a b x0 y0 α % --------------------------------- ellipse= [ 2 .92 .69 0 0 90; -.98 .8740 .6624 0 -.0184 90; -.02 .3100 .1100 .22 0 72; -.02 .4100 .1600 -.22 0 108; .01 .2500 .2100 0 .35 90; .01 .0460 .0460 0 .1 0; .01 .0460 .0460 0 -.1 0; .01 .0460 .0230 -.08 -.605 0; .01 .0230 .0230 0 -.606 0; .01 .0460 .0230 .06 -.605 90 ]; for k = 1:10 asq = ellipse(k,2)^2; % a^2 bsq = ellipse(k,3)^2; % b^2 phi = ellipse(k,6)*pi/180; % α,旋转角度 x0 = ellipse(k,4); % x 坐标 y0 = ellipse(k,5); % y 坐标 A = ellipse(k,1); % 折射指数 x=xg-x0; % 使椭圆居中 y=rot90(xg)-y0; % 数组旋转90度 cosp = cos(phi); sinp = sin(phi); idx=find(((x.*cosp + y.*sinp).^2)./asq + ((y.*cosp - x.*sinp).^2)./bsq |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |