0001
0002
0003
0004
0005
0006
0007
0008
0009 close all
0010 clear all
0011
0012 addpath(genpath('../General_Functions/'));
0013 addpath('Functions')
0014 addpath('Data')
0015
0016
0017
0018
0019 init_rand = 6;
0020 type_data = 10;
0021
0022
0023 resolution = 40;
0024
0025
0026
0027 type_robust = 2;
0028 out_C = 0;
0029
0030
0031 print_type = 0;
0032 plot_type = 0;
0033
0034
0035 ss = plot_init;
0036
0037
0038 init_rand_seed(init_rand);
0039
0040
0041
0042 display('----- Fig. 16.23: Reconstruction of the surface of a facade - dem from bundle adjustment result -----')
0043
0044
0045 [points,BB,delta_x,sigma_k,tmp] = ...
0046 simulate_points_dem_10('fa2_aurelo_result_pyra0_ausgeschnitten-1.ply',resolution);
0047
0048
0049 starttime = cputime;
0050 out_in = zeros(size(points,1),1);
0051 [ds,S,Sigma,Np,Nr,Mc,ver,A,w,w_f,W] = smooth_dem_robust_bilinear...
0052 (points,BB,delta_x,sigma_k,out_C,type_robust,type_data,out_in,print_type,plot_type);
0053 display([' complete time for solution: ',num2str(cputime-starttime)])
0054
0055
0056
0057 figure('name','Fig 16.23 dem from bundle adjustment result','color','w',...
0058 'Position',[0.5*ss(1),0.3*ss(2),0.4*ss(1),0.5*ss(2)]);
0059 plot_surface(ds,BB,delta_x,'plotfun',@mesh,'view',[12,30]);
0060 axis equal;axis off;
0061 title('fitted dem - min curvature^2','FontSize',16)
0062
0063