人脸识别之light_cnn

light_cnn出自2016 cvpr吴翔A Light CNN for Deep Face Representation with Noisy Labels，

优势在于一个很小的模型和一个非常不错的识别率。主要原因在于，

（1）作者使用maxout作为激活函数，实现了对噪声的过滤和对有用信号的保留，从而产生更好的特征图MFM(Max-Feature-Map)。这个思想非常不错，本人将此思想用在center_loss中，实现了大概0.5%的性能提升，同时，这个maxout也就是所谓的slice+eltwise，这2个层的好处就是，一，不会产生训练的参数，二，基本很少耗时，给人的感觉就是不做白不做，性能还有提升。

（2）作者使用了NIN(Network inNetwork)来减少参数，并提升效果，作者提供的A模型是没有NIN操作的，B模型是有NIN操作的，2个模型的训练数据集都是CASIA，但是性能有0.5%的提升，当然代价是会有额外参数的产生。但是相比其他网络结构，使用NIN还是会使模型小不少，作者论文中的网络结构和B,C模型相对应。

模型比较如下：

作者论文中的识别率是98.80%和实际测试结果非常接近，本人测试使用的对齐方式是使用作者提供的对齐方式。

         这里的S模型，是大神分享的基于light_cnn稍作修改的模型，其实就是通过卷基层替换全连接层来减少参数，当然为了保证输出为256维特征，作者在第一个卷基层的stride*2，这样使得卷积后特征图缩小为原来1半。这种方法减少了模型参数，同时，stride*2使得卷积过程少了1半以上，也许GPU并行后未必会有明显的反应，但是cpu模式下，就会看出运行时间减少不少。

 

修改的作者的对齐程序：

基本思想还是反射变换warpaffine的思想，根据眼睛做旋转，根据眼睛中点和嘴巴中点做缩放，然后crop而出）。

function [res, eyec2, cropped, resize_scale] = align(img, f5pt, ec_mc_y, ec_y,crop_size)f5pt = double(f5pt);ang_tan = (f5pt(1,2)-f5pt(2,2))/(f5pt(1,1)-f5pt(2,1));ang = atan(ang_tan) / pi * 180;img_rot = imrotate(img, ang, 'bicubic');imgh = size(img,1);imgw = size(img,2);% eye centerx = (f5pt(1,1)+f5pt(2,1))/2;y = (f5pt(1,2)+f5pt(2,2))/2;% x = ffp(1);% y = ffp(2);ang = -ang/180*pi;%{x0 = x - imgw/2;y0 = y - imgh/2;xx = x0*cos(ang) - y0*sin(ang) + size(img_rot,2)/2;yy = x0*sin(ang) + y0*cos(ang) + size(img_rot,1)/2;%}[xx, yy] = transform(x, y, ang, size(img), size(img_rot));eyec = round([xx yy]);x = (f5pt(4,1)+f5pt(5,1))/2;y = (f5pt(4,2)+f5pt(5,2))/2;[xx, yy] = transform(x, y, ang, size(img), size(img_rot));mouthc = round([xx yy]);resize_scale = ec_mc_y/abs(mouthc(2)-eyec(2));img_resize = imresize(img_rot, resize_scale);res = img_resize;eyec2 = (eyec - [size(img_rot,2)/2 size(img_rot,1)/2]) * resize_scale + [size(img_resize,2)/2 size(img_resize,1)/2];eyec2 = round(eyec2);img_crop = zeros(crop_size, crop_size, size(img_rot,3));% crop_y = eyec2(2) -floor(crop_size*1/3);crop_y = eyec2(2) - ec_y;crop_y_end = crop_y + crop_size - 1;crop_x = eyec2(1)-floor(crop_size/2);crop_x_end = crop_x + crop_size - 1;box = guard([crop_x crop_x_end crop_y crop_y_end], size(img_resize,1));if (box(2)>size(img_resize,2)||box(4)>size(img_resize,1))    img_crop(box(3)-crop_y+1:box(4)-crop_y+1, box(1)-crop_x+1:box(2)-crop_x+1,:) =imresize(img_resize,[box(4)-box(3)+1,box(2)-box(1)+1]);else    img_crop(box(3)-crop_y+1:box(4)-crop_y+1, box(1)-crop_x+1:box(2)-crop_x+1,:) = img_resize(box(3):box(4),box(1):box(2),:);end% img_crop = img_rot(crop_y:crop_y+img_size-1,crop_x:crop_x+img_size-1);cropped = img_crop/255;endfunction r = guard(x, N)x(x<1)=1;x(x>N)=N;r = x;endfunction [xx, yy] = transform(x, y, ang, s0, s1)% x,y position% ang angle% s0 size of original image% s1 size of target imagex0 = x - s0(2)/2;y0 = y - s0(1)/2;xx = x0*cos(ang) - y0*sin(ang) + s1(2)/2;yy = x0*sin(ang) + y0*cos(ang) + s1(1)/2;end

测试程序：

其中，list.txt数据格式如下：

路径/图片类别eye_x  eye_y  eye_x  eye_y  nose_x nose_y  mouse_x  mouse _y

clear;clc;imagelist=importdata('list.txt');addpath('/home/caffe/matlab');caffe.reset_all();% load face model and creat networkcaffe.set_device(0);caffe.set_mode_gpu();%-------1model = 'LightenedCNN_A_deploy.prototxt';weights = 'LightenedCNN_A.caffemodel';%-------2% model = 'LightenedCNN_B_deploy.prototxt';% weights = 'LightenedCNN_B.caffemodel';%-------3% model = 'LightenedCNN_C_deploy.prototxt';% weights = 'LightenedCNN_C.caffemodel';net = caffe.Net(model, weights, 'test');features=[];for k=1:size(imagelist.textdata,1)        f5pt = [imagelist.data(k,3),imagelist.data(k,5), imagelist.data(k,7),imagelist.data(k,9),imagelist.data(k,11);...        imagelist.data(k,2),imagelist.data(k,4), imagelist.data(k,6),imagelist.data(k,8),imagelist.data(k,10)];    f5pt=f5pt';       crop_size=128;    ec_mc_y=48;    ec_y=40;    [img2, eyec, img_cropped, resize_scale] = align(img, f5pt, ec_mc_y, ec_y, crop_size);        img_final = imresize(img_cropped, [crop_size crop_size], 'Method', 'bicubic');    img_final = permute(img_final, [2,1,3]);    img_final = img_final(:,:,[3,2,1]);        if size(img_final,3)>1    img_final = rgb2gray(img_final);    end        tic    res = net.forward({img_final});      toc    features=[features;[res{1,1}]'];   endscore=[];for i=1:2:size(features,1)    scoretmp=pdist2(features(i,:),features(i+1,:),'cosine');%chisq,emd,L1,cosine,euclidean    score=[score;scoretmp];endfigure,plot(score);caffe.reset_all();

训练的proto文件：http://download.csdn.net/detail/qq_14845119/9790137

reference：

https://github.com/AlfredXiangWu/face_verification_experiment