%% 数据标定，文件读取
clc,clear;
d = 0.2778;
nhalf = 512/2;
xc = -335/36; yc = 50/9; % 旋转中心在托盘的相应位置
a = xlsread('A题附件.xls','附件2','A1:FX512');
cc = xlsread('A题附件.xls','附件1','A1:IV256');
%% 添加高斯白噪声
a=awgn(a,0.005,5);%5表示噪声强度
%% 去噪模块
a=double(a);
h0=fspecial('gaussian',[5,5],1)
a=conv2(a,h0,'same');
%% CT反演模块
a = [zeros(100,180); a; zeros(100,180)];
img = iradon(a,[0:179]+120,'spline');
img=img*2.013;
n = size(img,1);%img矩阵的行数
[x, y] = meshgrid([-n/2:n/2]*d);%100/256
xp = [-50 50 50 -50 -50]';
yp = [-50 -50 50 50 -50]';
[xp, yp] = rotxy(xp, yp, -yc,-xc, 0);
imc=img(39:398,93:451);
imc=imresize(imc,[256,256]);
zj=ones(256,256);
for i=1:256
zj(:,256-i+1)=imc(i,:);%将图像顺时针旋转90度
end
jd=corr2(cc,zj)%原图像与反演图像相似度对比
