全卷積網(wǎng)絡(luò)FCN進行圖像分割

CNN能夠?qū)D片進行分類，可是怎么樣才能識別圖片中特定部分的物體，在2015年之前還是一個世界難題。神經(jīng)網(wǎng)絡(luò)大神Jonathan Long發(fā)表了《Fully ConvoluTIonal Networks for SemanTIc SegmentaTIon》在圖像語義分割挖了一個坑，于是無窮無盡的人往坑里面跳。

全卷積網(wǎng)絡(luò) Fully ConvoluTIonal Networks

CNN 與 FCN

通常CNN網(wǎng)絡(luò)在卷積層之后會接上若干個全連接層, 將卷積層產(chǎn)生的特征圖(feature map)映射成一個固定長度的特征向量。以AlexNet為代表的經(jīng)典CNN結(jié)構(gòu)適合于圖像級的分類和回歸任務(wù)，因為它們最后都期望得到整個輸入圖像的一個數(shù)值描述（概率），比如AlexNet的ImageNet模型輸出一個1000維的向量表示輸入圖像屬于每一類的概率(softmax歸一化)。

栗子：下圖中的貓, 輸入AlexNet, 得到一個長為1000的輸出向量, 表示輸入圖像屬于每一類的概率, 其中在“tabby cat”這一類統(tǒng)計概率最高。

FCN對圖像進行像素級的分類，從而解決了語義級別的圖像分割（semantic segmentation）問題。與經(jīng)典的CNN在卷積層之后使用全連接層得到固定長度的特征向量進行分類（全聯(lián)接層＋softmax輸出）不同，FCN可以接受任意尺寸的輸入圖像，采用反卷積層對最后一個卷積層的feature map進行上采樣, 使它恢復(fù)到輸入圖像相同的尺寸，從而可以對每個像素都產(chǎn)生了一個預(yù)測, 同時保留了原始輸入圖像中的空間信息, 最后在上采樣的特征圖上進行逐像素分類。

最后逐個像素計算softmax分類的損失, 相當(dāng)于每一個像素對應(yīng)一個訓(xùn)練樣本。下圖是Longjon用于語義分割所采用的全卷積網(wǎng)絡(luò)(FCN)的結(jié)構(gòu)示意圖：

簡單的來說，F(xiàn)CN與CNN的區(qū)域在把于CNN最后的全連接層換成卷積層，輸出的是一張已經(jīng)Label好的圖片。

其實，CNN的強大之處在于它的多層結(jié)構(gòu)能自動學(xué)習(xí)特征，并且可以學(xué)習(xí)到多個層次的特征：較淺的卷積層感知域較小，學(xué)習(xí)到一些局部區(qū)域的特征；較深的卷積層具有較大的感知域，能夠?qū)W習(xí)到更加抽象一些的特征。這些抽象特征對物體的大小、位置和方向等敏感性更低，從而有助于識別性能的提高。下圖CNN分類網(wǎng)絡(luò)的示意圖：

這些抽象的特征對分類很有幫助，可以很好地判斷出一幅圖像中包含什么類別的物體，但是因為丟失了一些物體的細節(jié)，不能很好地給出物體的具體輪廓、指出每個像素具體屬于哪個物體，因此做到精確的分割就很有難度。

傳統(tǒng)的基于CNN的分割方法：為了對一個像素分類，使用該像素周圍的一個圖像塊作為CNN的輸入用于訓(xùn)練和預(yù)測。這種方法有幾個缺點：一是存儲開銷很大。例如對每個像素使用的圖像塊的大小為15x15，然后不斷滑動窗口，每次滑動的窗口給CNN進行判別分類，因此則所需的存儲空間根據(jù)滑動窗口的次數(shù)和大小急劇上升。二是計算效率低下。相鄰的像素塊基本上是重復(fù)的，針對每個像素塊逐個計算卷積，這種計算也有很大程度上的重復(fù)。三是像素塊大小的限制了感知區(qū)域的大小。通常像素塊的大小比整幅圖像的大小小很多，只能提取一些局部的特征，從而導(dǎo)致分類的性能受到限制。

而全卷積網(wǎng)絡(luò)(FCN)則是從抽象的特征中恢復(fù)出每個像素所屬的類別。即從圖像級別的分類進一步延伸到像素級別的分類。

全連接層 -> 成卷積層

全連接層和卷積層之間唯一的不同就是卷積層中的神經(jīng)元只與輸入數(shù)據(jù)中的一個局部區(qū)域連接，并且在卷積列中的神經(jīng)元共享參數(shù)。然而在兩類層中，神經(jīng)元都是計算點積，所以它們的函數(shù)形式是一樣的。因此，將此兩者相互轉(zhuǎn)化是可能的：

對于任一個卷積層，都存在一個能實現(xiàn)和它一樣的前向傳播函數(shù)的全連接層。權(quán)重矩陣是一個巨大的矩陣，除了某些特定塊，其余部分都是零。而在其中大部分塊中，元素都是相等的。

相反，任何全連接層都可以被轉(zhuǎn)化為卷積層。比如，一個 K=4096 的全連接層，輸入數(shù)據(jù)體的尺寸是 7?7?512，這個全連接層可以被等效地看做一個 F=7,P=0,S=1,K=4096 的卷積層。換句話說，就是將濾波器的尺寸設(shè)置為和輸入數(shù)據(jù)體的尺寸一致了。因為只有一個單獨的深度列覆蓋并滑過輸入數(shù)據(jù)體，所以輸出將變成 1?1?4096，這個結(jié)果就和使用初始的那個全連接層一樣了。

全連接層轉(zhuǎn)化為卷積層：在兩種變換中，將全連接層轉(zhuǎn)化為卷積層在實際運用中更加有用。假設(shè)一個卷積神經(jīng)網(wǎng)絡(luò)的輸入是 224x224x3 的圖像，一系列的卷積層和下采樣層將圖像數(shù)據(jù)變?yōu)槌叽鐬?7x7x512 的激活數(shù)據(jù)體。AlexNet使用了兩個尺寸為4096的全連接層，最后一個有1000個神經(jīng)元的全連接層用于計算分類評分。我們可以將這3個全連接層中的任意一個轉(zhuǎn)化為卷積層：

針對第一個連接區(qū)域是[7x7x512]的全連接層，令其濾波器尺寸為F=7，這樣輸出數(shù)據(jù)體就為[1x1x4096]了。

針對第二個全連接層，令其濾波器尺寸為F=1，這樣輸出數(shù)據(jù)體為[1x1x4096]。

對最后一個全連接層也做類似的，令其F=1，最終輸出為[1x1x1000]

實際操作中，每次這樣的變換都需要把全連接層的權(quán)重W重塑成卷積層的濾波器。那么這樣的轉(zhuǎn)化有什么作用呢？它在下面的情況下可以更高效：讓卷積網(wǎng)絡(luò)在一張更大的輸入圖片上滑動，得到多個輸出，這樣的轉(zhuǎn)化可以讓我們在單個向前傳播的過程中完成上述的操作。

舉個栗子：如果我們想讓224×224尺寸的浮窗，以步長為32在384×384的圖片上滑動，把每個經(jīng)停的位置都帶入卷積網(wǎng)絡(luò)，最后得到6×6個位置的類別得分。上述的把全連接層轉(zhuǎn)換成卷積層的做法會更簡便。如果224×224的輸入圖片經(jīng)過卷積層和下采樣層之后得到了[7x7x512]的數(shù)組，那么，384×384的大圖片直接經(jīng)過同樣的卷積層和下采樣層之后會得到[12x12x512]的數(shù)組。然后再經(jīng)過上面由3個全連接層轉(zhuǎn)化得到的3個卷積層，最終得到[6x6x1000]的輸出((12 – 7)/1 + 1 = 6)。這個結(jié)果正是浮窗在原圖經(jīng)停的6×6個位置的得分！

面對384×384的圖像，讓（含全連接層）的初始卷積神經(jīng)網(wǎng)絡(luò)以32像素的步長獨立對圖像中的224×224塊進行多次評價，其效果和使用把全連接層變換為卷積層后的卷積神經(jīng)網(wǎng)絡(luò)進行一次前向傳播是一樣的。
Evaluating the original ConvNet (with FC layers) independently across 224x224 crops of the 384x384 image in strides of 32 pixels gives an identical result to forwarding the converted ConvNet one time.

如下圖所示，F(xiàn)CN將傳統(tǒng)CNN中的全連接層轉(zhuǎn)化成卷積層，對應(yīng)CNN網(wǎng)絡(luò)FCN把最后三層全連接層轉(zhuǎn)換成為三層卷積層。在傳統(tǒng)的CNN結(jié)構(gòu)中，前5層是卷積層，第6層和第7層分別是一個長度為4096的一維向量，第8層是長度為1000的一維向量，分別對應(yīng)1000個不同類別的概率。FCN將這3層表示為卷積層，卷積核的大小 (通道數(shù)，寬，高) 分別為 (4096,1,1)、(4096,1,1)、(1000,1,1)?？瓷先?shù)字上并沒有什么差別，但是卷積跟全連接是不一樣的概念和計算過程，使用的是之前CNN已經(jīng)訓(xùn)練好的權(quán)值和偏置，但是不一樣的在于權(quán)值和偏置是有自己的范圍，屬于自己的一個卷積核。因此FCN網(wǎng)絡(luò)中所有的層都是卷積層，故稱為全卷積網(wǎng)絡(luò)。

下圖是一個全卷積層，與上圖不一樣的是圖像對應(yīng)的大小下標，CNN中輸入的圖像大小是同意固定resize成 227x227 大小的圖像，第一層pooling后為55x55，第二層pooling后圖像大小為27x27，第五層pooling后的圖像大小為13*13。而FCN輸入的圖像是H*W大小，第一層pooling后變?yōu)樵瓐D大小的1/4，第二層變?yōu)樵瓐D大小的1/8，第五層變?yōu)樵瓐D大小的1/16，第八層變?yōu)樵瓐D大小的1/32（勘誤：其實真正代碼當(dāng)中第一層是1/2，以此類推）。

經(jīng)過多次卷積和pooling以后，得到的圖像越來越小，分辨率越來越低。其中圖像到 H/32?W/32 的時候圖片是最小的一層時，所產(chǎn)生圖叫做heatmap熱圖，熱圖就是我們最重要的高維特征圖，得到高維特征的heatmap之后就是最重要的一步也是最后的一步對原圖像進行upsampling，把圖像進行放大、放大、放大，到原圖像的大小。

最后的輸出是1000張heatmap經(jīng)過upsampling變?yōu)樵瓐D大小的圖片，為了對每個像素進行分類預(yù)測label成最后已經(jīng)進行語義分割的圖像，這里有一個小trick，就是最后通過逐個像素地求其在1000張圖像該像素位置的最大數(shù)值描述（概率）作為該像素的分類。因此產(chǎn)生了一張已經(jīng)分類好的圖片，如下圖右側(cè)有狗狗和貓貓的圖。

upsampling

相較于使用被轉(zhuǎn)化前的原始卷積神經(jīng)網(wǎng)絡(luò)對所有36個位置進行迭代計算，使用轉(zhuǎn)化后的卷積神經(jīng)網(wǎng)絡(luò)進行一次前向傳播計算要高效得多，因為36次計算都在共享計算資源。這一技巧在實踐中經(jīng)常使用，一次來獲得更好的結(jié)果。比如，通常將一張圖像尺寸變得更大，然后使用變換后的卷積神經(jīng)網(wǎng)絡(luò)來對空間上很多不同位置進行評價得到分類評分，然后在求這些分值的平均值。

最后，如果我們想用步長小于32的浮窗怎么辦？用多次的向前傳播就可以解決。比如我們想用步長為16的浮窗。那么先使用原圖在轉(zhuǎn)化后的卷積網(wǎng)絡(luò)執(zhí)行向前傳播，然后分別沿寬度，沿高度，最后同時沿寬度和高度，把原始圖片分別平移16個像素，然后把這些平移之后的圖分別帶入卷積網(wǎng)絡(luò)。

如下圖所示，當(dāng)圖片在網(wǎng)絡(luò)中經(jīng)過處理后變成越小的圖片，其特征也越明顯，就像圖像中顏色所示，當(dāng)然啦，最后一層的圖片不再是一個1個像素的圖片，而是原圖像 H/32xW/32 大小的圖，這里為了簡化而畫成一個像素而已。

如下圖所示，對原圖像進行卷積conv1、pool1后原圖像縮小為1/2；之后對圖像進行第二次conv2、pool2后圖像縮小為1/4；接著繼續(xù)對圖像進行第三次卷積操作conv3、pool3縮小為原圖像的1/8，此時保留pool3的featureMap；接著繼續(xù)對圖像進行第四次卷積操作conv4、pool4，縮小為原圖像的1/16，保留pool4的featureMap；最后對圖像進行第五次卷積操作conv5、pool5，縮小為原圖像的1/32，然后把原來CNN操作中的全連接變成卷積操作conv6、conv7，圖像的featureMap數(shù)量改變但是圖像大小依然為原圖的1/32，此時圖像不再叫featureMap而是叫heatMap。

現(xiàn)在我們有1/32尺寸的heatMap，1/16尺寸的featureMap和1/8尺寸的featureMap，1/32尺寸的heatMap進行upsampling操作之后，因為這樣的操作還原的圖片僅僅是conv5中的卷積核中的特征，限于精度問題不能夠很好地還原圖像當(dāng)中的特征，因此在這里向前迭代。把conv4中的卷積核對上一次upsampling之后的圖進行反卷積補充細節(jié)（相當(dāng)于一個差值過程），最后把conv3中的卷積核對剛才upsampling之后的圖像進行再次反卷積補充細節(jié)，最后就完成了整個圖像的還原。

缺點

在這里我們要注意的是FCN的缺點：
1、是得到的結(jié)果還是不夠精細。進行8倍上采樣雖然比32倍的效果好了很多，但是上采樣的結(jié)果還是比較模糊和平滑，對圖像中的細節(jié)不敏感。

2、是對各個像素進行分類，沒有充分考慮像素與像素之間的關(guān)系。忽略了在通常的基于像素分類的分割方法中使用的空間規(guī)整（spatial regularization）步驟，缺乏空間一致性。

實踐

輸入的圖片是：

現(xiàn)在可以直接來點代碼嗎？
# import package
import numpy as np
from PIL import Image
import caffe

# 初始化地址
caffe_root = fcn.berkeleyvision.org-master/voc-fcn8s/
model_def = caffe_root + deploy.prototxt # 模型文件
model_weights = caffe_root + fcn8s-heavy-pascal.caffemodel #模型權(quán)重值
test_image = images/2007_000129.jpg #測試圖片

# load image, switch to BGR, subtract mean, and make dims C x H x W for Caffe
im = Image.open(test_image)
in_ = np.array(im, dtype=np.float32)
in_ = in_[:,:,::-1] # change RGB image to BGR image
in_ -= np.array((104.00698793,116.66876762,122.67891434))
in_ = in_.transpose((2,0,1)) # Reshape the image from (500, 334, 3) to (3, 500, 334)

net = caffe.Net(model_def, model_weights, caffe.TEST) #導(dǎo)入模型
net.blobs[data].reshape(1, *in_.shape)
net.blobs[data].data[...] = in_ #讀入圖像
net.forward() #圖片進入前饋卷積神經(jīng)網(wǎng)絡(luò)
out = net.blobs[score].data[0].argmax(axis=0) #最后得到的圖片

print net.blobs[score].data[0].shape #(21, 500, 334)
print net.blobs[score].data[0].argmax(axis=0)

好了，已經(jīng)用fcn模型訓(xùn)練網(wǎng)一張圖片了，接著就是要看看圖片到底是怎么樣的楽
import matplotlib.pyplot as plt
# display plots in this notebook
%matplotlib inline

# set display defaults

print out.shape
plt.imshow(out)

輸出是：

現(xiàn)在做圖片分割的都是基于FCN的升級版、FCN超級升級版，F(xiàn)CN改版、FCN超級改版。。。個人覺得最難的、也是個人正在學(xué)習(xí)的是從如何研究自己的樣本，什么樣的樣本集才能提高最后結(jié)果的精度和召回率；有了樣本然后怎么給CNN訓(xùn)練，訓(xùn)練后如何把CNN->FCN，然后到FCN能夠?qū)π碌臄?shù)據(jù)進行分割。整套流程能夠自動化下來就更加perfect了。