Christensen G E, Joshi S C, Miller M I
Department of Electrical and Computer Engineering, The University of Iowa, Iowa City 52242, USA.
IEEE Trans Med Imaging. 1997 Dec;16(6):864-77. doi: 10.1109/42.650882.
This paper presents diffeomorphic transformations of three-dimensional (3-D) anatomical image data of the macaque occipital lobe and whole brain cryosection imagery and of deep brain structures in human brains as imaged via magnetic resonance imagery. These transformations are generated in a hierarchical manner, accommodating both global and local anatomical detail. The initial low-dimensional registration is accomplished by constraining the transformation to be in a low-dimensional basis. The basis is defined by the Green's function of the elasticity operator placed at predefined locations in the anatomy and the eigenfunctions of the elasticity operator. The high-dimensional large deformations are vector fields generated via the mismatch between the template and target-image volumes constrained to be the solution of a Navier-Stokes fluid model. As part of this procedure, the Jacobian of the transformation is tracked, insuring the generation of diffeomorphisms. It is shown that transformations constrained by quadratic regularization methods such as the Laplacian, biharmonic, and linear elasticity models, do not ensure that the transformation maintains topology and, therefore, must only be used for coarse global registration.
本文展示了猕猴枕叶和全脑冷冻切片图像的三维(3-D)解剖图像数据以及通过磁共振成像获得的人类大脑深部脑结构的微分同胚变换。这些变换以分层方式生成,兼顾了全局和局部解剖细节。初始的低维配准是通过将变换限制在低维基中来完成的。该基由放置在解剖结构中预定义位置的弹性算子的格林函数和弹性算子的本征函数定义。高维大变形是通过模板与目标图像体积之间的不匹配生成的矢量场,该矢量场被约束为纳维 - 斯托克斯流体模型的解。作为该过程的一部分,跟踪变换的雅可比行列式,以确保生成微分同胚。结果表明,受诸如拉普拉斯、双调和及线性弹性模型等二次正则化方法约束的变换不能确保变换保持拓扑结构,因此,只能用于粗略的全局配准。