【数学建模】CUMCM-2017A CT系统参数标定及成像 思路及部分代码

前面

蛮遗憾的一次建模……怎么说呢,在当时的能力下可以做的更好的。
题附下面git库了……。感觉前两问回答的比较正常第4问相当于没怎么写……

思路

第四问太糟糕了不管了……

1

第一题的成像图
这里写图片描述
然后把两条曲线分离…没找到好的方法直接用用图像处理……
这里写图片描述

然后可以根据这两个情况算出每组数据对应的角度(这里应该可以用到数值分析的方法……我是直接大步长粗搜再小步长……就是那么暴力……
拟合之后得到每次的旋转角度。

2&3

2和3的思路差不多……
然后先把第一题的数据代进去验证正确性然后再扔2&3的数据(这样会好一点
滤波反投影其实比反投影重建算法就多了个滤波。
写迭代重建算法是我傻……(还是最后一个晚上写的,没学过线代的人一边写一边MDZZ 结果发现并不好用……
然后我们居然还撒fufu的分析了一下各种不同滤波器的效果【队友闲着没事干
这里写图片描述

这里写图片描述

这里写图片描述

然后这里最大的问题是有些点会出现负数看着就很难过orz

当时处理是,对于第二问,把轮廓找出来了,然后再对里面的图案弄了个分水岭算法把它分块(因为第二问明显是比较规则的轮廓
取轮廓(其实这个也写了很久……将近6h……被队友嫌弃死……
这里写图片描述
分水岭之后分了区域的图
这里写图片描述

Code

都是matlab码
像我这种代码习惯差的半年之后看到自己命名奇怪的码内心是炸的……
只记得

function [img,H] = relvbo(varargin)

这个是滤波反投影的函数orz

https://github.com/birdy-C/CUMCM2017A.git

后面

刚刚在找相关资料的时候看到的……

如何评价2017数学建模国赛? - 鸽子的迷信的回答 - 知乎
https://www.zhihu.com/question/65346481/answer/248032685

对了参考文献的话
这里写图片描述
【1】这个是很基础的的,我记得提到把离散回到连续的方法
【3】主要对着这个写的 是一本 很破 很旧 年纪比我大的 还附有借书卡的书!(深深感觉到自己的姿势水平落后了时代十几年

猜你喜欢

转载自blog.csdn.net/birdy_/article/details/79405770