Herzberg William, Rowe Daniel B, Hauptmann Andreas, Hamilton Sarah J
Department of Mathematical and Statistical Sciences; Marquette University, Milwaukee, WI 53233 USA.
Research Unit of Mathematical Sciences; University of Oulu, Oulu, Finland and with the Department of Computer Science; University College London, London, United Kingdom.
IEEE Trans Comput Imaging. 2021;7:1341-1353. doi: 10.1109/tci.2021.3132190. Epub 2021 Dec 2.
The majority of model-based learned image reconstruction methods in medical imaging have been limited to uniform domains, such as pixelated images. If the underlying model is solved on nonuniform meshes, arising from a finite element method typical for nonlinear inverse problems, interpolation and embeddings are needed. To overcome this, we present a flexible framework to extend model-based learning directly to nonuniform meshes, by interpreting the mesh as a graph and formulating our network architectures using graph convolutional neural networks. This gives rise to the proposed iterative Graph Convolutional Newton-type Method (GCNM), which includes the forward model in the solution of the inverse problem, while all updates are directly computed by the network on the problem specific mesh. We present results for Electrical Impedance Tomography, a severely ill-posed nonlinear inverse problem that is frequently solved via optimization-based methods, where the forward problem is solved by finite element methods. Results for absolute EIT imaging are compared to standard iterative methods as well as a graph residual network. We show that the GCNM has good generalizability to different domain shapes and meshes, out of distribution data as well as experimental data, from purely simulated training data and without transfer training.
医学成像中大多数基于模型的学习图像重建方法都局限于均匀域,比如像素化图像。如果基础模型是在非均匀网格上求解,这源于非线性逆问题典型的有限元方法,那么就需要进行插值和嵌入操作。为克服这一问题,我们提出了一个灵活的框架,通过将网格解释为图并使用图卷积神经网络来制定我们的网络架构,从而将基于模型的学习直接扩展到非均匀网格。这就产生了所提出的迭代图卷积牛顿型方法(GCNM),该方法在逆问题的求解中包含了正向模型,而所有更新都是由网络直接在特定问题的网格上计算得出的。我们展示了电阻抗断层成像的结果,这是一个严重不适定的非线性逆问题,通常通过基于优化的方法来求解,其中正向问题通过有限元方法求解。将绝对电阻抗断层成像的结果与标准迭代方法以及图残差网络进行了比较。我们表明,GCNM对不同的域形状和网格、分布外数据以及实验数据具有良好的通用性,这些数据来自纯模拟训练数据且无需迁移训练。