目录
一、任务分析
二、题目一求解
题目分析：
模型建立：
模型求解：
3.1读数据，绘图
3.2探测器之间的距离计算
3.3小圆与发射—接收系统的相对运动
3.4提取小圆圆心位置坐标
三、题目二求解
一、任务分析
历年赛题
任务：CT标定，CT图像重建
知识：图像处理
标定模板为椭圆与球。CT测量所得数据为512*180的二维矩阵。
行数512:接收点个数，列数180:CT旋转次数。所成的二维矩阵并不是真实物体的成像，而是不断旋转的水平投影的叠加。
每个角度都有一个大小为512的投影向量，一共180个角度的叠加。
二、题目一求解
求解：CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离、CT系统使用的X射线的180个方向
题目分析：
标定 CT 系统的安装误差
模型建立：
（1）探测器之间的距离
模型求解：
3.1读数据，绘图
表格一数据读取
A=xlsread('E:\math_model_data\data\A.xls');imshow(A)
表格二数据读取
B=xlsread('E:\math_model_data\data\A.xls',2);X=1:size(B,2);Y=1:size(B,1);figure(1)mesh(X,Y,B)
3.2探测器之间的距离计算
B=xlsread('E:\math_model_data\data\A.xls',2);figurelenb = (1:length(B));subplot(221)scatter(lenb( B(:,1)~=0 ),B( B(:,1) ~= 0 ,1),20,'filled')title('1')subplot(222)scatter(lenb( B(:,6)~=0 ),B( B(:,6) ~= 0 ,6),20,'filled')title('6')subplot(223)scatter(lenb( B(:,11)~=0 ),B( B(:,11) ~= 0 ,11),20,'filled')title('11')subplot(224)scatter(lenb( B(:,60)~=0 ),B( B(:,60) ~= 0 ,60),20,'filled')title('60重合')
图1:球投影29个点。图6:球投影29个点。 图11:球投影28个点。
解得，探测器之间的距离和小圆直径的不等关系式:
由图即可解得近似解：
3.3小圆与发射—接收系统的相对运动
模板：小球与椭圆，小球投影为小圆。系统：CT系统。
模板与系统相对运动，一动一静。以系统为参照物，系统看做静止物体，则模板为运动物体。取第一个探测器为坐标原点，建立坐标系，解算
根据特殊点解得、与
点A:
点B:
点C：
解得
3.4提取小圆圆心位置坐标
分离段：取极值
重合段：边缘提取后，取平均值作为中心
得到180个小圆圆心，解得180个角度。
3.5参数修正
3.5.1 进一步确定旋转角度的范围
3.5.2 修正旋转中心坐标
三、题目二求解
模型求解：
3.1输出重建图像
滤波反投影重建算法
