(17)图像分割——U-Net

    U-Net: Convolutional Networks for Biomedical  Image Segmentation

    PDF :https://arxiv.org/pdf/1505.04597.pdf

(1)网络结构

    U-Net是一个生物图像分割的网络,因为其形状像U型,得名U-Net。


U-Net网络结构(最少32*32的分辨率)

    U-net网络非常简单,前半部分作用是特征提取,后半部分是上采样。这样的结构也叫做编码器-解码器结构。U-net与其他常见的分割网络有一点非常不同的地方:U-net采用了完全不同的特征融合方式:拼接,U-net采用将特征在channel维度拼接在一起,形成更厚的特征。而FCN融合时使用的对应点相加,并不形成更厚的特征。

    已有的研究中如Ciresan,取以图像中每个像素点为中心的局部区域(patch),进行训练,相对于单张图片增大了数据集。但是存在以下的缺点:1、由于每个patch之间可能存在很多重叠,造成了重复训练;2、由于全连接层的使用,大的patch需要更多的最大池化层,而最大池化层丢失信息严重,特别是位置准确率。因此考虑不同尺度特征图的信息可以同时兼顾位置和内容。

    U-Net左半部分采用2个3*3的卷积(无填充)接一个Relu和一个2*2的最大池化(步长为2),每次下采样后,通道数翻倍。右半部分通过上采样增大特征图的分辨率,为了精准定位,左右两边对应大小的特征图进行拼接,增大了特征图的通道,可以使内容信息传播到更深的层中(这里主要数考虑细胞边缘等信息,之前的下采样造成了很多信息损失,在高分辨率和更抽象特征当中做一个折衷),以获得更准确的分割。

    最后用过1*1 的卷积将通道数压缩到分类的数目。整个网络结构中一共23个卷积层。由于没有全连接层,因此可以输入任意尺度(满足下采样要求的)的图像进行端到端的训练。这里示意图中最终分为2类,例如第一张表示的是第一类的得分(即每个像素点对应第一类都有一个得分),第二张表示第二类的得分heatmap,然后作为softmax函数的输入,算出概率比较大的softmax类,选择它作为输入给交叉熵进行反向传播训练

(2)损失函数:

    能量函数(the energy function)是通过与交叉熵损失函数相结合的最终特征图,并利用像素级的soft-max函数来计算的。加权交叉熵定义为:

    E= \sum_{x\in \Omega }w(x)log(p_{l(x)} (x))其中,l={1,...,K}表示每个像素正确的标签,权重函数w。这里相当于最大似然估计,求最大值,这里的损失函数是求最大值,也就是无限接近于0的值,因为p(x)值是大于0小于1的,这里的log肯定就是小于0 的,所以0就是它的最大值。

    soft-max函数定义为:


a k (x)表示在第x像素点上第k个特征通道的激活函数。K为类别数,p k (x)为近视最大函数(the approximated maximum-function)。例如Pk(x) 约等于1表示ak(x)取得最大激活,而其他通道上的Pk(x)则近似为0

    权重函数定义为:


其中,Wc 表示平衡类别频率的权重图,d1表示到最近细胞边界的距离,d2表示到第二近细胞边界的距离。基于经验,将w 0 设为10,σ约为5像素。 即在细胞边界附近的像素点给的权重会大一些,离细胞比较远的像素点的权重会小一些。因为,如果同类细胞贴的比较近,可能就会增大训练的难度,减少准确率,毕竟卷积会考虑该像素点周围的一些特征,而两个相同的类的细胞贴在一起,就容易误判,所以对这种两个相同类贴在一起的细胞边界,给予较大的权重,使的训练之后分类分割更准确。


    权重初始化基于标准差为√(2/N)的高斯分布,其中N表示每个神经元的输入结点数。例如3*3的卷积核,64通道时N = 3*3*64 = 576 

    训练过程使用带动量的随机梯度下降,动量大小为0.99。

(3)数据扩增:

        在小的数据集上做了大量的数据扩增,由于数据量比较小,我们使用了弹性形变的方法进行数据扩展,对于医疗数据来说变形是生物组织最常见的变异,通过这种方法可以对变形进行模拟。这个是3*3的网格上使用随机位移矢量产生平滑形变,其中位移来自于10像素标准差的高斯分布,且通过双三次插值法计算得出。具体操作原理见弹性形变的博客。这里的label通过形态学的方法获得,看做原图的一个通道,同时做形变等操作。

    因此U-Net适用于小规模的数据集,同时取得比较准确的分割。而U-Net也适合于医疗数据。因为医疗数据和其他数据不同,没有用ImageNet进行预训练。对于边界上的像素信息,残缺信息通过镜像手段进行补全。

(4)实验和效果


图像分割效果

    通过对照实验,并在多个数据集上训练,取得了最好的效果。这里为了最大效率利用GPU,减少了patch的大小,同时增加了每个输入图像的大小,这里一个patch就是一个tile(输入图像):注:在切面图像中,一整个切片是(slice),一个slice里剪切出多个tile。

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