qgao21 / CoreDiff

35 stars 3 forks source link

关于系统矩阵 #8

Open mkxyz-123 opened 7 months ago

mkxyz-123 commented 7 months ago

你好!我在阅读有关迭代算法的文献的时候看到了系统矩阵这个概念,成像的过程就可以写为AX=P(其中,A为系统矩阵;X为图像数据;P为得到的投影数据;在这里X为一维的列向量,得到的P也是一维的列向量)。如果要使用梯度下降法来重建图像的话,那么目标函数就可以写成||AX-P||^2。我在编写相关的代码,因此需要得到系统矩阵A。我知道A中的元素 aij 代表图像中第 j 个像素对第 i 个投影值的贡献,但是我没有找到具体得出系统矩阵A的算法。我在网上搜到的帖子基本都是获得投影数据P或者反投影的方法,没有对这个中间过程A的描述。请问你有没有接触过这些?如果你知道有什么算法可以得出A可以告诉我吗(或者一些相关的知识)。非常感谢!!!

qgao21 commented 7 months ago

系统矩阵A的计算方式有四种:像素驱动、 射线驱动、距离驱动以及面积驱动,其中像素驱动和射线驱动方法应用较多。具体计算公式相关文献有提及,你也可以通过MIRT工具包来构建生成系统矩阵。

你留个联系方式吧,微信或者邮箱都可以,后面交流方便。

mkxyz-123 commented 7 months ago

好的!感谢!我的微信号是shangsiluting