平均人脸的制作探讨

smiles54

已知：有10幅已截好的人脸bmp图片，种族肤色各不相同，图片规格都是48*48
问：求10个人的平均脸并显示出来，如何编程？
我的想法：
for i=1:10
      a=imread(strcat('E:\平均脸\pjl',num2str(i),'.bmp'));
      b=a(1:48*48);
% b是行矢量 1×N，其中N＝48*48，提取顺序是先列后行，即从上到下，从左到右
      b=double(b);
      allsamples=[allsamples; b];  
% allsamples 是一个M * N 矩阵，allsamples 中每一行数据代表一张图片，M为10
end
samplemean=mean(allsamples); % 平均图片，1 × N
imshow(samplemean);
问题：显示不出图像，说是数据太大

eight

原帖由 smiles54 于 2007-9-13 00:10 发表

                               
登录/注册后可看大图

已知：有10幅已截好的人脸bmp图片，种族肤色各不相同，图片规格都是48*48
问：求10个人的平均脸并显示出来，如何编程？
我的想法：
for i=1:10
      a=imread(strcat('E:\平均脸\pjl',num2str(i),'.bmp')) ...

使用 imread 后，a 矩阵就是二维（灰度图）或者三维（RGB图）的，不明白你为何要做成一维的

stephenhope

恩 首先你要确定你样本是不是都是灰度图。
还有你的 samplemean=mean(allsamples); % 平均图片，1 × N
要是你要imshow的话 也要resize 成 48*48才行啊！

smiles54

谢谢楼上的看法，以上是我最初的设想，是把十幅图都读到一个矩阵中，再求平均，现在我直接读入10幅图像，再求平均即可显示成了，但是图像只是数据的平均，并不是所有图像特征的融合。
现在我用wavemenu里的fuse工具已经可以实现两个图像的融合了，但是十幅图像还不能实现，如果用for循环，量会有影响，不知该怎么对图像的质着手。以下是平均法和融合法分别得到的人脸，看看效果先。

[ 本帖最后由 smiles54 于 2007-9-14 22:08 编辑 ]