addpath('toolbox_signal') addpath('toolbox_general') addpath('toolbox_graph') addpath('solutions/meshproc_3_denoising') clear options; name = 'nefertiti'; name = 'elephant-50kv'; options.name = name; % useful for displaying [X0,F] = read_mesh(name); n = size(X0,2); m = size(F,2); options.lighting = 1; % clf; plot_mesh(X0,F, options); rho = 0.015; N = compute_normal(X0,F); X = X0 + repmat(rho*randn(1,n),[3,1]).*N; clf; plot_mesh(X,F,options); axis('tight'); E = [F([1 2],:) F([2 3],:) F([3 1],:)]; E = unique_rows([E E(2:-1:1,:)]')'; E0 = E(:,E(1,:)