环形伪影的添加

转载 2018年04月15日 13:55:07

基于Shepp-Logan添加环形伪影
clc
clear all
L = 256;%图像大小
P = phantom(L);
theta = 0:359;
R = radon(P,theta);%获得正弦图

for ii = 120:120
    for jj = 1:360
        R(ii,jj) = 0.0;
    end
end
%figure,imshow(R,[])
I = iradon(R, theta, 'linear', 'Hamming', L);%%获得重建图像
JG = 2*pi/(4*L);%%角度采样频率
width = 2*pi/(JG) + 1;
[m,n] = size(I);
radius = ceil(sqrt(m^2 + n^2)/2);
POL = zeros(radius*2, width);
for gama = 0:JG:(2*pi)
    for rad = 0.25:0.25:radius   %%半径采样频率
        i = round(gama/JG+1);%%round四舍五入取整
        j = (rad-0.25)*4 + 1;
        x = rad*cos(gama);%%x坐标变换
        y = rad*sin(gama);%%y坐标变换
        ceil_x = ceil(x);%ceil向下取整
        ceil_y = ceil(y);
        floor_x = floor(x);
        floor_y = floor(y);
        chae3 = abs(x - floor_x);
        chae4 = abs(y - floor_y);
        x1 = floor_x + L/2;
        x2 = ceil_x + L/2;
        y1 = floor_y + L/2;
        y2 = ceil_y + L/2;
        if x1>L
            x1 = L;
        end
        if x1<1
            x1 = 1;
        end
        if x2>L
            x2 = L;
        end
        if x2<1
            x2 = 1;
        end
        if y1>L
            y1 = L;
        end
        if y1<1
            y1 = 1;
        end
        if y2>L
            y2 = L;
        end
        if y2<1
            y2 = 1;
        end
        POL(j,i) = (1-chae3)*(1-chae4)*I(x1,y1) + chae3*(1-chae4)*I(x2,y1) + (1-chae3)*chae4*I(x1,y2) + chae3*chae4*I(x2,y2);%%双线性插值
    end
end
figure, hold on
subplot(1,3,1),imshow(I)%直接重建图像
subplot(1,3,2),imshow(POL)%极坐标变化后图像
subplot(1,3,3),imshow(R)%正弦图

Unity 3D残影技术实现

主要实现了3D游戏中残影效果,非常适合游戏开发。
  • 2015年05月23日 15:48

CT重建图像环状伪影的去除方法

  • 2017年08月09日 23:26
  • 1.17MB
  • 下载

sobel算子边缘检测去除环形伪影

  • 2015年05月19日 21:59
  • 1KB
  • 下载

shepp-logan模型笛卡尔坐标变换到极坐标,环状伪影仿真

利用shepp-logan模型进行环状伪影仿真,处理方法为后处理方法,即对CT图像进行处理,将图像从笛卡尔坐标变换到极坐标系中。clc clear all L = 256;**图像大小 P = pha...
  • wangkun7673
  • wangkun7673
  • 2017-08-13 20:23:19
  • 399

canny算子边缘检测去除环形伪影

  • 2015年05月19日 21:51
  • 1021B
  • 下载

CT图像金属伪影校正算法研究

  • 2011年05月26日 10:07
  • 4.41MB
  • 下载

【Cocos2d-x】自定义环形菜单

思路 1.定义一个类继承CCMenu,并重写itemForTouch方法,该方法在CCMenu被触摸的时候调用,用于获取触摸的菜单条目。在itemForTouch方法中自定义触摸判定。 2.定义...
  • linchaolong
  • linchaolong
  • 2015-03-08 20:57:42
  • 1379

软件工程(十三)

过程设计语言(PDL)PDL也称为伪码。 如:if I>0 then 执行订单数据输入模块 els...
  • u013553804
  • u013553804
  • 2016-05-18 00:44:29
  • 408

显微CT常用名词解释

原文http://www.xjgkyy.com/gksy/show.asp?id=657 CT值    CT值(CT number)是以水的CT值为零,而相对于其他物质X线的衰减值。例如,空气的...
  • MNaivety
  • MNaivety
  • 2011-04-22 15:01:00
  • 1782

MFC的sheep logan图像的CT扫描重建

  • 2010年03月29日 19:43
  • 362KB
  • 下载
收藏助手
不良信息举报
您举报文章:环形伪影的添加
举报原因:
原因补充:

(最多只允许输入30个字)