林叶很快就进入了状态。</p>
吃饱了饭，又喝了一杯咖啡提神，此刻精神十分亢奋。</p>
“离散模型中使用代数迭代法（ART）对网格化的离散数据直接建模，用矩阵迭代的方法对像素网格进行处理，实现图像的重构。”</p>
这是林叶的思路。</p>
至于数学模型具体如何，林叶还在思考。</p>
这个离散模型，最重要的还是迭代问题。</p>
FBP（滤波反投影）算法具有的缺陷与限制，只能应用代数迭代才能完美解决。</p>
“这里得先引入投影矩阵的概念...。”</p>
林叶嘀咕说道。</p>
时间一点一点的过，林叶空白的A4纸上现在已经写满了数学公式。</p>
r_ij={1，i号射线通过j号像素中心//0，其他。</p>
记x=[x1，x2，...，xJ]^T为J维图像矢量，x_j表示图像上第j个像素点的吸收强度，</p>
p=[p1，p2...pI]^T为I维投影矢量，p_i表示第i条射线所经过的所有像素的总射线投影;</p>
Rx=p;</p>
...</p>
林叶看着草稿纸上最为核心的迭代模型，</p>
x^(k+1)=x^(k)+λ^(k)(q_ikr_ikx^(k)/||r_ik||^2)</p>
林叶嘀咕道：</p>
“这里还差一点东西，那就是降噪算法。不论何种重建算法，不可避免会带入伪影点。</p>
组委会给的数据含有噪声，成像质量可能并不会达到最好的要求。这个其余学校的建模选手也应该会认识到这个问题才对。”</p>
“这里必须得用算法进行降噪。我得查一查这方面的文献。”</p>
林叶用笔压着笔记本电脑前的草稿纸，然后查阅起了文献。</p>
“这种问题应该在国外医学SCI期刊上比较好寻找合适的算法。”</p>