十、成本函数

为了估计转换模型的哪些参数可以最好地同时显示两个图像,我们需要一种方法来定义图像之间的差异,这被称为成本函数。当图像对准良好时,一个好的成本函数应该很小,并且随着它们逐渐变得越来越不对准,一个好的成本函数应该变得更大。选择合适的成本函数在很大程度上取决于您尝试配准的图像类型。如果图像是相同类型的(例如,跨不同时间点重新对准fMRI数据),则成本函数只需要确定跨两个图像的图像强度的相似性。如果图像完全对齐,则图像中的强度值应该彼此非常接近(暂时不考虑它们可能由于激活等有趣因素而改变的事实)。这个问题通常被称为“医疗模式内”配准。另一方面,如果图像具有不同类型的对比度(例如,T1加权的MRI图像和T2加权的图像),则最优对齐将不会在图像之间产生相似的值。这被称为“模态之间”配准。对于T1加权图像和T2加权图像,在T1加权图像中,白质将比灰质亮,而在T2加权图像中,灰质比白质亮(参见图2.4),因此我们不能简单地匹配所有图像的强度。相反,我们希望使用一种对不同体素集的相对强度敏感的方法;例如,我们可能希望将一幅图像的亮部与另一幅图像的暗部进行匹配。

2.4 不同MRI图像类型的示例。不同脑区(如白质、灰质、脑室)的相对强度因图像类型的不同而不同,这意味着它们不能简单地通过匹配图像间的强度来对齐。

在这里,我们描述了几个最常见的成本函数,用于MRI图像的模态内和模态间配准。

1.  最小平方

最小二乘成本函数可能是最熟悉的,因为它是大多数标准统计方法的基础。此成本函数衡量每个图像中体素强度之间的均方差:

其中A_{v} B_{v} 分别指图像A和B中的第v个体素的强度。因为它测量每个体素上的值的相似性,所以最小二乘代价函数仅适用于模态内配准。即使在模态中,如果两幅图像具有不同的强度分布(例如,其中一幅整体比另一幅更亮或强度范围更广),它的性能也会很差。一种方法是在使用最小二乘成本函数之前首先对强度分布进行缩放,这是AIR软件包中的一个选项,以便它们在图像上落在相同的范围内。

2. 归一化相关

归一化相关性测量两幅图像中体素强度之间的线性关系。它被定义为:

这一措施仅适用于医疗机构内的配准。在进行运动校正的许多不同成本函数的比较中,詹金森(Jenkinson,2002年)发现,归一化相关性比包括最小二乘在内的其他几种成本函数的配准更准确。它是FSL软件包中运动校正的默认成本函数。

2.3 互信息

尽管前面描述的用于模态内配准的成本函数以经典统计学为基础,但互信息成本函数(Pluim等人,2003)(可用于模态间或模态内配准)源自信息论中的熵概念。熵指的是信号中存在的不确定性或随机性的量:

其中p_{i} 是变量的每个可能值x_{i} 的概率;对于连续变量,这些值被分组到N个箱中,通常称为直方图箱。熵测量变量的每个不同可能值在信号中出现的程度。如果只能出现一个信号值(即,xi的概率为pi=1),熵将被最小化。如果每个不同的值出现的频率相等(即,对于所有Xi,pi=1/N),则熵最大化。以这种方式,它与信号的方差以及不确定性密切相关,人们可以利用该不确定性来预测信号的下一个值。通过检查图像的联合直方图,可以将熵扩展到多个图像,该直方图绘制出图像中所有体素的所有可能值上强度组合的频率(参见图2.5)。如果两个图像相同,则联合直方图仅在对角线上具有值(因为在每个图像中的体素s中的值将是相同的),而图像之间的差异导致值在直方图上的更大分散;请注意,对于这种模态内的情况,相关性将是比互信息(MI)更合适的成本函数度量。对于来自不同模式的图像,在互信息更合适的情况下,配准错误越大,联合直方图中的离散度越大(见图2.6)。然后,可以从该联合直方图计算两个图像A和B的联合熵,如下所示:

其中i为A的索引值,j为B的索引值,pij为两者同时出现的概率。当图像B的值完全可由图像A中的相同体素的值预测时,该度量最低。

互信息是各个图像的总熵和联合熵之间的差:

其中H(A)和H(B)是分别为每个图像中的值单独计算的熵(称为边际熵),H(A,B)是联合熵。当联合熵最小时,互信息最大,当一幅图像的值可从另一幅图像中最大限度地预测时就会发生这种情况。因此,互信息可以用作两个图像之间的相似性的度量。

互信息的一个潜在问题是,在某些情况下,即使图像之间的重叠减少,互信息也会增加。为此,建议对互信息系数进行归一化(Studholme et al., 1999):

所有主要的软件包(FSL、SPM和AFNI)都提供了用于图像配准的常规和归一化互信息代价函数。

2.5 T1加权图像和T2加权图像之间的联合直方图示例。图像的暗度表明直方图的这一区间出现的频率更高。显示了每幅图像中两个不同体素的强度,一个在灰质(顶部切片),另一个在白质(底部切片)。

2.6 图2.5中T1加权和T2加权图像的联合直方图。左面板显示了原始配准图像的联合直方图。中间面板显示其中一个图像旋转1度后的联合直方图,而右侧面板显示其中一个图像旋转180度后的联合直方图。每个比较的MI值都显示在面板标签中;MI随着图像变得不那么相似而减小。

2.4 相关比

相关性比率(Roche et al.,1998)衡量一个度量中的方差被另一个度量中的方差捕获的程度。两个图像A和B的相关比定义为

其中,k是B的每个唯一值的索引,N是B的唯一值的数目。如果A和B相同,则在具有某些特定B值的所有体素上,A值将没有方差,并且相关比变为零。该方法类似于最初在AIR软件包中为PET-MRI共同配准实施的Woods标准(Woods等人,1993年),尽管它在某些情况下表现不同(Jenkinson&Smith,2001)。它既适用于模态内配准,也适用于模态间配准,是FSL软件包中默认的模态间成本函数。


以上内容来自《Handbook of functional MRI Data Analysis》。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,732评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,496评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,264评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,807评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,806评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,675评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,029评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,683评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 41,704评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,666评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,773评论 1 332
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,413评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,016评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,978评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,204评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,083评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,503评论 2 343

推荐阅读更多精彩内容