(基于粘滞流体B样条模型的快速非刚体配准方法.docVIP

基于粘滞流体B样条模型的快速非刚体配准方法 荣成城1) 周健2) 曹国刚3) 罗立民4) 1) 3) 4)（东南大学影像科学与技术实验室, 中国南京, 210096.） 2)（法国国家健康研究院，信号与图像处理实验室，法国雷恩，35000.）分别为待配准的源图像和目标图像函数，其中表示像素空间坐标。如前文所述，配准的目的是计算如下的位移场。 我们可以定义如下的相似能量泛函： 然后通过最小化该能量（1）来获得相应所需的位移场。在粘性流体配准方法中，u是通过求解偏微分方程来得到的。让我们来考虑单个流体粒子的运动过程。设位于处的粒子时刻所受外力为，它是关于当前位置和位移的函数。设该粒子受力而产生的瞬时速度为，则根据粘滞流体运动学原理[9]，则有如下偏微分方程： 此处为拉普拉斯算子，为梯度算子，为散度运算，为粘滞系数，本文考虑一般情况，取，。该偏微分方程反映了u随v和k变化而变化的关系。此时，k可以取泛函关于u的变分形式： 因此，最优的位移场u可以通过求解偏微分方程（2）来获得。 通过引入时域的离散化操作，第一个偏微分方程可以改写为： 其中表示时间间隔。 在二维空间中展开和，可以得到： 通过离散化v1，v2，k1，k2，Christensen【4】等使用逐次超松弛来求解上述方程组。 2.2引入B样条 求解偏微分方程组（5）的准确性取决于有限