2017数学建模总结 ) |4 Y) U$ I" w1 R4 K7 f' U# b. N5 D% \
迟到的建模总结,昨晚12点交完论文,打印材料,今天一天满课,刚刚把上周作业写完的我,终于有空闲时间写写这个总结了。 0 U% O5 E% q8 \% I 9 Z2 o. h u+ ~: \5 x, X$ K3 A简单概括游戏体验非常差!!! 1 K) X. m. p: \2 t7 n, h+ f0 A+ D# O. L% g: {6 t0 f' S7 `
第一天晚上选题,看了看A,看了看B,B题分析性太强,属于应用统计,数据分析类,和队友商量着战略性放弃,然后就开始查资料,搞东西了,突然回忆起以前姜起源的书上貌似讲过CT成像这种东西,翻出来看果然有这玩意儿。然后看了看题目发现精度是4位小数,似乎可以暴力枚举答案啊(ACMer的特性),再仔细对照姜老的模型,发现参数太多了,180个角度,还有中心,连单元距离也没告诉,,,,图片是256x256,几乎搜到建模结束也不可能搜出来,突然陷入绝望,本着第一晚看看题的想法,就这样回宿舍碎觉了。) B, o' Q+ q; z. H" _/ y. k
. O H0 z( w$ ]% e0 j5 k第二天上午我突然发现了一个很好的几何特性,圆无论从那边投影,始终间距是不变的,那不就是数一数探测器格子的事吗,于是乎敲了几行代码,发现180次投影绝大多数都是29个格子,那就8/29了。然后发现由于椭圆中兴到圆心的距离知道,而中心峰值最大,那就那么通过这个就可以大致确定角度了,还想了想其实中兴也可以这样算啊,取探测器中间位置为中心线不就行了?几何关系就弄出来了,到底靠不靠谱啊,不敢想,这时大概10.30,唐老进来吹逼了,还以为它能告诉我们一个很好的思路,结果就是瞎逼逼,连CT成像的数字图像处理原理都不知道……radon变换都没提到,间距都不会求。哎感到无望的我决定先去把自己的想法实现.敲了一敲,角度还是大致精确,对着其实角度29.。。。几度一个一个往上跑,可是中兴当取不同连线的时候就不行了.又陷入绝望,2,3两问还要重建图像呢,到了下午也不知道怎么处理,这时一个老师道出了精髓,你怎么就能确定中心一定在256像素呢。瞬间茅塞顿开.再也找不到求中心的方法了,在晚上要结束的时候发现matlab有iradon变换的函数,随便放了108个角度进去,竟然出来图像了,虽然有点倾斜,然后搜了一搜,发现scipy里还是准备好这些图像处理的基本变换函数的,就接着有python试了试,果然也是可以的., x( x% k2 K# ]- A! o1 k* G! D