十二、再细分和插值

一旦估计了用于变换的参数,然后将它们应用于原始图像以创建变换后的或重新分割的图像。这包括使用变换后的图像坐标以及原始强度填充坐标空间中每个体素的值。如果变换受到限制,使得变换的位置与原始图像中的位置完全重叠,则可以简单地使用来自相应变换的体素的值填充这些体素。然而,一般而言,变换将涉及体素的分数;例如,在运动校正中,被校正的运动通常小于体素尺寸的1/10。在这种情况下,变换后的体素与原始体素并不完全重叠,因此有必要对原始图像强度的值进行插值,以获得重新拼接的图像。

1. 最近邻插值

在最近邻插值中,新体素的值将替换为原始图像中最近体素的值。这种形式的interpolation很少使用,因为它存在许多问题。当对同一图像执行多次内插时,这一点尤其明显,如图2.9所示。其次,在使用最近邻内插时,变换参数的连续变化可能会导致成本函数值的不连续变化,这使得它不适合用于优化方法(通常假设成本函数是连续的)。

优化最近邻内插的一种情况是变换其中体素值表示标签而不是物理强度的图像。例如,某些软件包中包含的解剖图集被存储为具有特定但任意值的图像,这些值对应于图集中的每个结构(例如,海马体素的值为12,而杏仁核中的体素的值为20 )。使用在给定无意义结果的这些数字之间取平均值的插值方法,将使用邻域插值来确保变换图像中的标签值保持与原始图像中的值相同的精确值。

2.9 这是重复插值效果的一个极端例子。每幅图像以0.01弧度旋转六次,并在每一步分别使用最近邻、三线性和正弦插值进行重新分割。请注意,与sinc插值相比,最近邻图像中产生伪影和三线性图像的模糊度增加。

2. 线性插值

这种方法在应用于三维时通常被称为三线性插值,它涉及对原始图像中每个最近点处的值进行加权平均。图2.10显示了线性内插的一个例子。与高阶插值方法相比,该方法具有相对较快的优点,因为它只考虑紧邻新位置的那些点。但是,与诸如Sinc插值之类的高阶插值方法相比,它倾向于使图像更加模糊。

2.10 线性插值的一个例子。点a、b、c和d是原始图像位置,位于间距为2 mm的网格上。点g等效于沿X维度向左平移0.75 mm、沿Y维度向下平移0.75 mm后的点b。为了确定图像在点g处的强度,首先沿一维进行插值;在本例中,我们沿X维进行插值,找到点e和f的值。然后,我们沿着另一个轴进行插值,以确定点g处的值。

3. 高阶插值

与最近邻居(仅使用单个最近体素)和线性内插(在三维中集成八个最近体素)相比,已经开发了许多内插方法来集成跨越更广泛的体素集合的信息。最常见的高阶插值方法是sinc插值,它使用sinc函数[sinc(x)=sin(x)/x],如图2.11所示。原则上,这种形式的插值应该使用来自图像中每个体素的信息,因为sinc函数扩展到无穷大。这在计算上会非常昂贵。通过使用加窗SINC函数使SINC插值更可行,其中该函数仅从被插值点延伸有限距离,而不是覆盖整个图像。有许多不同形式的窗口可以应用;常见选项包括汉宁窗和矩形窗。汉宁窗口似乎导致相对于矩形窗口的内插误差减小(Ostuni等人,1997年),因此如果可用,应该使用半径(或半长)至少为四个体素的窗口来选择。

2.11 Sinc函数的一个示例。X轴上的刻度线表示原始栅格位置,插值位置处的刻度线为0。为了获得插值值,将在原始网格上的每个位置放置SINC函数(由图像强度缩放);然后通过将新位置处的所有SINC函数值相加来获得插值值。(这相当于使用sinc函数对图像进行卷积。)

另一种形式的高阶插值使用基函数,用于空间转换模型。基函数,提供了一种更一般的插值方法,它可以包括最近邻插值、线性插值以及高阶非线性插值。


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

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

推荐阅读更多精彩内容

  • 网上有很多介绍插值算法的,但感觉收获都不大 简书的编译器可能有点问题tex公式在预览模式公式可以显示,从主页浏览就...
    zht1999阅读 3,417评论 0 0
  • opencv interpolation 选项 所用的插值方法INTER_NEAREST 最近邻插值IN...
    reallocing阅读 898评论 0 1
  • 常见插值算法 在现实生活中,我们经常会遇到把图像进行放大、几何空间变换的情况等等,这些操作都需要在源图像和目标图像...
    瞎了吗阅读 5,343评论 0 1
  • 一、图像的像素级运算 点运算 点运算具有如下特点:1)点运算针对图像中的每一个像素灰度,独立地进行灰度值的...
    老子不吃鱼_阅读 3,814评论 0 0
  • 本章介绍了基于elastix的基本配准概念。 更高级的配准主题将在第6章中讨论。图像配准是医学影像领域的重要工具。...
    peterpan_hai阅读 9,677评论 1 10