4 L: [$ k( B: D9 f- `- K
目录
一、任务分析
二、题目一求解
题目分析：
模型建立：
模型求解：
3.1读数据，绘图
3.2探测器之间的距离计算
3.3小圆与发射—接收系统的相对运动
3.4提取小圆圆心位置坐标
三、题目二求解
B3 E! x# F0 z
一、任务分析历年赛题
任务：CT标定，CT图像重建
知识：图像处理
标定模板为椭圆与球。CT测量所得数据为512*180的二维矩阵。
行数512:接收点个数，列数180:CT旋转次数。所成的二维矩阵并不是真实物体的成像，而是不断旋转的水平投影的叠加。
每个角度都有一个大小为512的投影向量，一共180个角度的叠加。
二、题目一求解求解：CT系统旋转中心在正方形托盘中的位置
、探测器单元之间的距离
、CT系统使用的X射线的180个方向
题目分析：标定 CT 系统的安装误差
模型建立：（1）探测器之间的距离
模型求解：3.1读数据，绘图表格一数据读取
: k* l! K* w5 z1 x; F. l" s
4 b2 l X/ U0 o" q) PA=xlsread('E:\math_model_data\data2017\A.xls');
3 h; L' v' _' V$ ^8 b, R& O* D8 J3 W: U2 ]3 Z8 J7 V% Z" K; _
- 6 E' g5 H3 C; D) t0 x) a, D4 {
5 K y' B, A/ W# P! G3 I; n
imshow(A)
. q% _7 ?8 K3 J6 G1 R0 E* E- d
4 w: m4 a" |. L+ a7 C1 ^
4 f& x) k9 B+ @/ F3 z
3 v) W; U1 V( h
表格二数据读取
- + {" t' J; \, M& P0 E- W/ H- m
& ?) W5 i! h6 B3 z. XB=xlsread('E:\math_model_data\data2017\A.xls',2);
1 f1 } m6 q/ q; s. w- L- C# w" V+ K8 W4 t
- 4 [/ E# w& p2 |: w
$ Y1 A# @# f8 Y: NX=1:size(B,2);
4 n( f% w7 c$ O0 a! N/ \8 J$ u$ n2 ~: h
- 7 X: Q$ ~1 H( }6 m, ]
z3 P% y0 ]2 L4 I. sY=1:size(B,1);* ~, J# c# ~. q2 o7 @
. L5 \- v$ p* f, L$ X
( ?$ Y3 Y: b) E+ N/ K! T+ l5 u6 W2 b' F
figure(1)
4 v: ~) {, i% |4 i) L* x* g1 l. l6 E% y5 u) K3 Q9 [7 U2 ^
) A+ a! n" \; ?! i3 |8 i4 O7 k
mesh(X,Y,B)
- Q& m; D- V! V8 ^ h0 J
! Z3 c+ Z" W* x* g5 O; |$ ], i" ?$ H3 f5 ?! ~
- a6 W/ \) F8 N5 F# I
3.2探测器之间的距离计算
6 r3 F' @7 y0 p8 |! P, g7 @6 T3 C4 R7 \
B=xlsread('E:\math_model_data\data2017\A.xls',2);5 h& p' J1 C n( j6 E; V+ {' n
0 \& j- a: ?/ B& B* u' t! ]# ?' X
: c; p3 \ [" G: k9 Z+ j( z
7 C3 Q3 n9 U# X& [0 h
5 U, C5 p! @ E, U0 i9 |
" A- o& G1 K( F; \- 5 W6 w$ d& w8 H t. u( S% i8 g
0 ?" j/ Y: d' r& ~' R1 b* C8 M2 ufigure& M1 {! a2 k! W1 i
+ P z& M0 p' |% H; u R+ x - 6 M5 k. W9 f# K( W) z6 k( C% l
1 T" \/ `2 u5 i9 m8 Klenb = (1:length(B));
, @- { w4 r0 \" F* i5 S3 C& e2 H) G/ T% u# { D B3 t
# N/ y/ w* q, M6 R/ e1 s5 V. e
% C/ ?0 e) h' k P) X- n8 U, _. p. x, n7 ^
# s4 T, q6 b) \' J2 f
' {9 f' D8 |- S" _- |
- X, K- E9 v6 I/ isubplot(221)
% f- J- p+ x5 h7 R( u7 s" D, h/ ]( T& ^" m4 t2 E& U2 n& F5 ~ o
- 0 V( D; @* d9 ^/ H; N a t
1 d# E5 Y3 U% g; g
scatter(lenb( B(:,1)~=0 ),B( B(:,1) ~= 0 ,1),20,'filled')8 R0 t& d, [7 V4 B; Z- G# T2 ^
9 W% ^$ x0 q4 {! |7 v - - n2 A3 \6 Q; Z9 {3 b
) G# d, _* j3 J d+ U2 _( r) R! Q
title('1')
_6 C- ^" z4 }. t; q$ ` o. S* V8 p0 |+ {3 y& E
* e1 y' d/ `" h% |
8 H8 i( x: m$ v7 a$ D, N8 U
2 s5 \1 T! }' q6 [& _3 ]- ~- Z2 m; F
- 6 h2 y' J; [; f! B8 q& Q- W) c
+ n. y, v" v5 q, X4 rsubplot(222)' I2 @' G; ]3 e a0 x$ c1 M
, C# _' [4 h/ m V: u - ! N c3 m3 \5 |* W4 L
1 F P( S" w. }6 y O' t/ ^scatter(lenb( B(:,6)~=0 ),B( B(:,6) ~= 0 ,6),20,'filled')
. B& ]$ o, c# k# S- F( }' E! q
; o0 h T5 O% y4 w# A
$ r2 @6 Z* A" S) P" n# y3 l$ s- R9 V" h5 a
title('6') l+ Y% Z1 C5 B: j
" D- Q; T! s) a7 b, X$ q
; a* U" u& ~$ |0 a+ k6 C2 {
0 t0 J4 G4 c! l
9 I+ y T+ @9 S' A/ y7 c8 w( e& q7 o" P
- - F' \6 Q v' ~" j
" h u% }4 Y( q, A& J. o& Y8 q
subplot(223)
3 K2 H! T5 H! @3 B" c: o: k% S
% [3 ^- Y; o) \" }6 ?" N - 7 i" _9 n- C! _* o; L
, p$ n$ L9 O8 ~2 @7 R! bscatter(lenb( B(:,11)~=0 ),B( B(:,11) ~= 0 ,11),20,'filled')
) F& y$ Q- B0 N' a, M+ a, U; d
9 n0 Q+ a5 p) c- } - & s; I: H x4 i2 ]4 t1 r- c
\ S3 P1 t* S2 C4 K
title('11')
: T9 z' `5 @( `; d( o' x0 @
) @# T5 P5 Y) v8 U; g- [# z - 6 b4 N: `: j; s* e
4 s) w: O" h' Z$ B
' h P# V8 Z5 Z; q" `& v V) ?% G/ y* Y9 \
- 4 m) \; `& ~0 U5 s0 ?
4 a1 R: t) c6 O \8 i( E1 A$ Zsubplot(224)# E; S* ?1 e. S) }! `, w
' b. {7 s5 B/ K
- 7 M" z' _! p( X9 {4 Y
7 ], F0 R7 c, J+ w
scatter(lenb( B(:,60)~=0 ),B( B(:,60) ~= 0 ,60),20,'filled')
3 c/ m' C/ w+ G# j( `& S+ _2 e5 u3 c2 h
; N# C/ A! W8 ^5 l$ K0 D1 W! E' q
5 \3 f9 }% r* Q- W; a1 s3 Mtitle('60重合')
+ ^5 ]" y/ a, ?( ]! ?% }: P, U4 ~* V& G/ H! K+ ~
5 y/ }8 b" X. {2 k3 h
图1:球投影29个点。图6:球投影29个点。 图11:球投影28个点。
解得，探测器之间的距离和小圆直径的不等关系式: 
由图即可解得近似解：
3.3小圆与发射—接收系统的相对运动模板：小球与椭圆，小球投影为小圆。系统：CT系统。
模板与系统相对运动，一动一静。以系统为参照物，系统看做静止物体，则模板为运动物体。取第一个探测器为坐标原点，建立坐标系，解算
根据特殊点解得
、
与
; U: j2 ^8 X& y& E1 Z6 F1 _
点A:
点B:
点C：
解得
3.4提取小圆圆心位置坐标分离段：取极值
6 |! u; V% ?* x, |0 r4 }4 W. _1 ] r重合段：边缘提取后，取平均值作为中心
得到180个小圆圆心，解得180个角度。
3.5参数修正3.5.1 进一步确定旋转角度的范围
, n; z& M% J2 A7 M
3.5.2 修正旋转中心坐标