$ @, s- v2 b. E% r$ x; {第二天上午我突然发现了一个很好的几何特性,圆无论从那边投影,始终间距是不变的,那不就是数一数探测器格子的事吗,于是乎敲了几行代码,发现180次投影绝大多数都是29个格子,那就8/29了。然后发现由于椭圆中兴到圆心的距离知道,而中心峰值最大,那就那么通过这个就可以大致确定角度了,还想了想其实中兴也可以这样算啊,取探测器中间位置为中心线不就行了?几何关系就弄出来了,到底靠不靠谱啊,不敢想,这时大概10.30,唐老进来吹逼了,还以为它能告诉我们一个很好的思路,结果就是瞎逼逼,连CT成像的数字图像处理原理都不知道……radon变换都没提到,间距都不会求。哎感到无望的我决定先去把自己的想法实现.敲了一敲,角度还是大致精确,对着其实角度29.。。。几度一个一个往上跑,可是中兴当取不同连线的时候就不行了.又陷入绝望,2,3两问还要重建图像呢,到了下午也不知道怎么处理,这时一个老师道出了精髓,你怎么就能确定中心一定在256像素呢。瞬间茅塞顿开.再也找不到求中心的方法了,在晚上要结束的时候发现matlab有iradon变换的函数,随便放了108个角度进去,竟然出来图像了,虽然有点倾斜,然后搜了一搜,发现scipy里还是准备好这些图像处理的基本变换函数的,就接着有python试了试,果然也是可以的." E, x4 i% p& [3 E
2 l& w7 E# a/ Q7 C, |0 j5 B/ A& \第三天就是各种方法乱上了,因为iradon没有中心这个接口,自己又不会写,反变换像素与原图像素不一致,各种乱搞,发现还原出来的图居然正好将原图缩小的一些,对比探测器长度和和原图宽度,竟然是根号2倍关系,随机明白了,开始截图,还原像素了,就这样度过了惊人的第3天. 0 s4 A6 I: c9 ] d3 O; u+ [7 t, }第四天分析精度再也弄不出来什么东西了,只能xjb搞了,老师说几何方法在4问肯定没法做,那我也没办法了……可是老师说枚举180个角度确定中心去最小误差的搜索方法我至今怀疑,怎么可能一个角度就能搜出中心,而且,一次iradon变换运行时间至少为1秒啊,(可能是搞ACM的原因本能觉得是不可行的)。 8 Z. A0 R) @$ W6 P/ I! g5 t! C- i/ v- K/ P
在晚上的时候老唐居然在群里发出了,最恶心的嘲讽,附上截图 $ K8 J# Z' q9 K, n: F) H
) g; j2 Q+ k d. o6 J0 j' Q
# R( K8 i1 W& \5 Z
估计做A题的人的心得凉透了….就这样在这种嘲讽中忙活过了最后的4小时改论文时间…..8 R% a( J A+ N, |
2 Q, }# X8 S: ]7 E2 V# J
至此,数学建模经历全剧终。其实从暑假培训起就有点不想搞了,培训就xbb,什么也没讲到,原理从来不讲,就叫你用matlab命令搞一搞就行了,全程用python的我根本用不到这些东西的,而且现在python社区如此庞大,本科生搞搞小科研个人认为完全够用了,不知道学校为何还要让大家集体用盗版的matlab。(机房也是盗版的).记得当时听见有老师说建模是他们的事情,我们所要做的就是编程实现,,,,,,这次估计打脸了,在第二天,大佬唐,(就是那个最后挖苦我们的人)说我们知道的也就那么多了,其余我们也不会了….然后一直xbb。。# Z! \7 s6 |- n% X+ t6 s' _# U; b# P
--------------------- ' @/ s4 x7 O5 P: `5 ~
作者:taotao.zou 4 W$ q- n" J& b5 I, r4 w4 z, t9 y/ V
来源:CSDN e5 [* L" e! b' R* _; n + Q' T0 X8 `) L/ F* i/ B- G. E8 R% x; H3 @+ x8 m9 {) u+ y4 L