우리는받은 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)