opencv实现多角度重建

根据多角度拍摄图像进行三维重构 能得出三维模型

function ptch = build_model（ voxels ）
% build_model（ voxels ）
%  This function is used to draw a surface based on the provided voxels 
%  structure. It uses Matlab‘s ISOSURFACE command to do that.
% 
% ARGUMENTS:
% VOXELS = this argument contains the voxels to be drawn. 
%
% RETURNS:
% PTCH = SHOWSURFACE（VOXELS） also returns handles to the patches created.
%

% First puts the data in a grid
data_x = unique（voxels.XData）;
data_y = unique（voxels.YData）;
data_z = unique（voxels.ZData）;

% Then expands the model in each direction by one step
data_x = [data_x（1）-voxels.Resolution; data_x; data_x（end）+voxels.Resolution];
data_y = [data_y（1）-voxels.Resolution; data_y; data_y（end）+voxels.Resolution];
data_z = [data_z（1）-voxels.Resolution; data_z; data_z（end）+voxels.Resolution];

% Convert to a grid
[XYZ] = meshgrid（ data_x data_y data_z ）;

% Create a voxel grid with empty spaces and fill only those elements 
% where voxels are present.
V = zeros（ size（ X ） ）;
N = numel（ voxels.XData ）;
for ii=1:N
    ix = （data_x == voxels.XData（ii））;
    iy = （data_y == voxels.YData（ii））;
    iz = （data_z == voxels.ZData（ii））;
    V（iyixiz） = voxels.Value（ii）;
end

% Now draw it
ptch = patch（ isosurface（ X Y Z V 0.5 ） ）;
isonormals（ X Y Z V ptch ）
set（ ptch‘FaceColor‘ ‘black‘ ‘EdgeColor‘ ‘red‘ ）;

set（gca‘DataAspectRatio‘[1 1 1]）;
xlabel（‘X‘）;
ylabel（‘Y‘）;
zlabel（‘Z‘）;
view（130-25）
lighting（ ‘gouraud‘ ）
axis（ ‘tight‘ ）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-04-18 14:40  project\
     文件        1443  2013-05-31 19:17  project\build_model.m
     目录           0  2017-04-18 14:40  project\data\
     文件      388594  2013-05-31 19:17  project\data\0.jpg
     文件      395946  2013-05-31 19:17  project\data\1.jpg
     文件      496935  2013-05-31 19:17  project\data\10.jpg
     文件      488137  2013-05-31 19:17  project\data\11.jpg
     文件      441134  2013-05-31 19:17  project\data\12.jpg
     文件      427088  2013-05-31 19:17  project\data\13.jpg
     文件      458722  2013-05-31 19:17  project\data\14.jpg
     文件      525791  2013-05-31 19:17  project\data\15.jpg
     文件      582925  2013-05-31 19:17  project\data\16.jpg
     文件      502297  2013-05-31 19:17  project\data\17.jpg
     文件      511186  2013-05-31 19:17  project\data\18.jpg
     文件      548797  2013-05-31 19:17  project\data\19.jpg
     文件      361701  2013-05-31 19:17  project\data\2.jpg
     文件      548873  2013-05-31 19:17  project\data\20.jpg
     文件      518973  2013-05-31 19:17  project\data\21.jpg
     文件      550512  2013-05-31 19:17  project\data\22.jpg
     文件      597930  2013-05-31 19:17  project\data\23.jpg
     文件      368549  2013-05-31 19:17  project\data\3.jpg
     文件      315052  2013-05-31 19:17  project\data\4.jpg
     文件      337683  2013-05-31 19:17  project\data\5.jpg
     文件      358522  2013-05-31 19:17  project\data\6.jpg
     文件      421499  2013-05-31 19:17  project\data\7.jpg
     文件      521621  2013-05-31 19:17  project\data\8.jpg
     文件      507293  2013-05-31 19:17  project\data\9.jpg
     文件        2039  2013-05-31 19:17  project\data\morp.mat
     文件        1466  2013-05-31 19:17  project\data\morp2.mat
     文件        1073  2013-05-31 19:17  project\generate_silhouettes.m
     文件         912  2013-05-31 19:17  project\generate_voxels.m
............此处省略7个文件信息