十年磨一剑

the format of .sel file is like:# per vertex status: 0=fixed, 1=deformable-region, 2=handle# 5261 vertices1111211101111function [group0, group1, group2] =

先讲下各个公式公式(2)是laplacian coordinates的定义公式(5)第一部分是保证laplacian coordinate坐标的一致性, 后面是handle点的约束公式(8)是T的构造公式(12)是求T中的系数2D MATLAB 版在它的project里面自带, 这里将其拓展成3D版本,  源代码如下:func

比较low的方法for i = 1:nC(:,:,i) = A(:,:,i) * B(:,:,i);endelegant的方法, 用arrayfunarrayfun(@(M1,M2,n) M1(:,:,n)*M2(:,:,n), A, B, size(A,3) )

arrayfun(@(n) norm(A(n,:)), 1:size(A,1))注意当返回不是标量的时候记得将'UniformOutput'参数改为0arrayfun(@(n) norm(A(n,:)), 1:size(A,1), 'UniformOutput', 0);

c0 = arrayfun(@(x) mean(v0([x neighbors{x}],:)), 1:n, 'UniformOutput', 0); c0 = reshape(cell2mat(c0),3,[]);c0 = c0'

function c = scalar2color(scalar, precision)%precision表示有多少间隙, % convert scalar to color vector%%% Copyright (c) 2017 seamanj@NCCAif size(scalar,2)>size(scalar,1) scalar = scalar';endc

function faces = perform_faces_reorientation(vertex,faces, options)% perform_faces_reorientation - reorient the faces with respect to the center of the mesh%% faces = perform_faces_reorientation

function [points,D] = perform_farthest_point_sampling_mesh( vertex,faces, points, nbr_iter, options )% perform_farthest_point_sampling - samples points using farthest seeding strategy%% [points,D]

function options = my_plot_mesh(vertex,face,options)% plot_mesh - plot a 3D mesh.%% my_plot_mesh(vertex,face, options);%% 'options' is a structure that may contains:% - 'normal' : a (

function show_initilization%show_target(0.2, 1)%show_template('v4_mesh0.obj', 0.5, 1, 'interp')name = 'target.obj';%options = [];options.name = name;[vertex,faces] = read_mesh(name);vert

%# figurefigure, set(gcf, 'Color','white')Z = peaks; surf(Z); axis tightset(gca, 'nextplot','replacechildren', 'Visible','off');[az,el]=view;xl=xlim;yl=ylim;zl=zlim;%# preallocatenFrames =

for frame = 55:-1:0 before_name = strcat(strcat('Q:\CGI2017\bak\RenderByHoudini\v4_mesh',num2str(frame)),'.obj');% add photo prefix after_name = strcat(strcat('Q:\CGI2017\bak\RenderByHoudin

function M = mat3D2blk(A, d)M = reshape(permute(A,[1 3 2]), [],d);M = mat2cell(M,(size(M,1)/d)*ones(2,1),d);M = blkdiag(M{:});

function main pics = dir('coil-20-proc');%list all the files in folder photo name = extractfield(pics, 'name');% extract name fields into a cell array index = ~ismember(name,[ {'.'},{'

[vertices1, faces1] = read_mesh('dog6.obj');[vertices1, faces1] = check_face_vertex(vertices1, faces1);X = vertices1(1,:);color = (X - min(X))./ (max(X) - min(X));color = color';colormap j...