十年磨一剑

原创 differential forms

Tutorial comes from https://www.youtube.com/watch?v=jFe6dMnpQho 0-forms : f:R3→Rf : \mathbb R^3 \rightarrow \mathbb R 1-forms : α=adx+bdy+cdz=w<a,b,c> \alpha =adx + bdy + cdz = w_{<a,b,c>} 2-forms

原创 compute_voronoi_triangulation_mesh

function [vertex_voronoi,faces_voronoi] = compute_voronoi_triangulation_mesh(Q, vertex, faces)% compute_voronoi_triangulation_mesh - compute a triangulation%% [vertex_voronoi,face_voronoi] = com

转载 perform_faces_reorientation

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

转载 perform_farthest_point_sampling_mesh

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]

原创 AE的一些基本知识

ctrl+shift+d切割帧composition setting 设置帧数add to render queue,  output to render 渲染视频到本地文件

原创 export frames as pictures

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

原创 write frames to a video

%# 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 =

转载 how to fix the frame size?

The problem is that due to the rotation of the graph the size of the axis changes, so you end up with a sequence of frames of different sizes. You can see this if you insert the liness = size(mov(k)

原创 sort_nat

function [cs,index] = sort_nat(c,mode)%sort_nat: Natural order sort of cell array of strings.% usage: [S,INDEX] = sort_nat(C)%% where,% C is a cell array (vector) of strings to be sorted.%

原创 how change files in matlab2 , cell 函数的运用

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,[ {'.'},{'

原创 intersectLineMesh3d

C:\TestMatlab\geom3dfunction [points, pos, faceInds] = intersectLineMesh3d(line, vertices, faces, varargin)%INTERSECTLINEMESH3D Intersection points of a 3D line with a mesh%% INTERS = intersectLi

原创 Deformable 3D shape registration based on local similarity transforms

本质上是shape matching的as similar as possible版本它的主要贡献在于, 找到一种方法使displacement field尽量平滑, 具体的做法就是公式7, 使每个点的偏移量, 尽量接近它周围邻居偏移量的平均值.注意tkt_k公式在paper中是错误的 应该为 tk=ck(t)−skRkC0kt_k = c_k(t) - s_kR_kC^0_k因为rk=skRk

原创 求cluster的质心坐标

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

原创 Laplacian surface editing

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

原创 houdini 渲染wireframe

原创 文章索引加1

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

原创 图片转视频软件

VirtualDub

原创 my plot mesh

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 (

原创 sub figure without letter numbering

\begin{figure*} \captionsetup[subfigure]{labelformat=empty} \centering \subfloat[Asdaf]{\includegraphics[width=0.11\textwidth]{./human2gorilla/crop/sample.png}\label{fig:sample}}\hspace{0.05\te

原创 read tetgen file

function [p,f,t] = read_tetgenmesh( file )% READ_TETGENMESH Reads Tetgen meshes% Reads the native Tetgen format mesh-files into Matlab. The output% [p,f,t] (nodes,faces,tetras) is similar to

原创 compute_vertex_ring

function vring = compute_vertex_ring(face)% compute_vertex_ring - compute the 1 ring of each vertex in a triangulation.%% vring = compute_vertex_ring(face);%% vring{i} is the set of vertices

原创 houdini 做选点效果

注意 object merge的时候, 选好transform

原创 如何算三角形的cotangent

公式为:cotC=a2+b2−c24A\cot C = \frac{a^2 + b^2 - c^2}{4A} 证明 根据余弦定理cosC=a2+b2−c22ab\cos C = \frac{a^2 + b^2 - c^2}{2ab} 根据面积公式 A=12absinCA = \frac{1}{2}ab \sin C 所以 sinC=2Aab\sin C = \frac{2A}{ab}所以co

原创 两个3*3*n旋转矩阵在第三维相乘

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

原创 scalar2color

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

原创 rescale

function y = rescale(x,a,b)% rescale - rescale data in [a,b]%% y = rescale(x,a,b);%% Copyright (c) 2004 Gabriel Peyr�if nargin<2 a = 0;endif nargin<3 b = 1;endm = min(x(:));M

转载 normrow

function [ B ] = normrow( A ) % NORMROW Compute l2 row vector norms % % B = normrow( A ) % % Input: % A #A by D list of row vectors of dimension D % Output: % B #A list of norms o