人脸识别之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 center
- x = (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;
- end
- function r = guard(x, N)
- x(x<1)=1;
- x(x>N)=N;
- r = x;
- end
- function [xx, yy] = transform(x, y, ang, s0, s1)
- % x,y position
- % ang angle
- % s0 size of original image
- % s1 size of target image
- x0 = 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 network
- caffe.set_device(0);
- caffe.set_mode_gpu();
- %-------1
- model = '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}]'];
- end
- score=[];
- for i=1:2:size(features,1)
- scoretmp=pdist2(features(i,:),features(i+1,:),'cosine');%chisq,emd,L1,cosine,euclidean
- score=[score;scoretmp];
- end
- figure,plot(score);
- caffe.reset_all();
训练的proto文件:http://download.****.net/detail/qq_14845119/9790137
reference:
https://github.com/AlfredXiangWu/face_verification_experiment