天天看點

matlab中文本框的邊框旋轉,ParaView Clip/Slice,Box參數:位置,旋轉,縮放的含義;在Matlab中重建框...

我們複制框屬性出paraview包的:

box = struct();

box.translate = [-274.8975, -114.0316, -333.6671];

box.rotate = [27.066, 119.62, -175.472];

box.scale = [1.031, 0.663, 1.2233];

箱尺寸是相對于剪輯過濾器上定義的幾何尺寸。是以,首先擷取幾何的範圍,其頂點存儲在surfpts中。

surfpts = ...; % [nPoints x 3]

mn = min(surfpts,[],1);

mx = max(surfpts,[],1);

在這個例子中,我想填充一個網格點的框。

box.xyz = arrayfun(@(a,b) linspace(a,b,n)', mn, mx, 'UniformOutput',false);

box.xyz = cell2mat(box.xyz);

[X,Y,Z] = ndgrid(box.xyz(:,1), box.xyz(:,2), box.xyz(:,3));

box.XYZ = [X(:), Y(:), Z(:)];

現在我們要根據我們做的盒子的Paraview 3D插件的動作來縮放,旋轉和平移我們的網格(順序!)。

% Scaling

box.XYZ = bsxfun(@times, box.XYZ, box.scale);

% Rotation

box.XYZ = (rotz(box.rotate(3)) * rotx(box.rotate(1)) * roty(box.rotate(2)) * box.XYZ')';

% Translation

box.XYZ = bsxfun(@plus, box.XYZ, box.translate);

最後,我們通過在ParaView中再次打開建立的框來檢查結果。它與盒子完全比對。

TR = delaunayTriangulation(box.XYZ);

[box_tri, box_pts] = TR.freeBoundary();

vtkwrite('tmp_box.vtk','polydata','triangle',box_pts(:,1),box_pts(:,2),box_pts(:,3),box_tri)