2017-04-19 19 views
0

3D 상자을 지정된 삼각형 모양의 지오메트리 옆에 특정 크기와 위치로 정의하고자합니다.ParaView 클립/슬라이스, 박스 매개 변수 : 위치, 회전, 축척의 의미. Matlab의 재구성 상자

따라서 ParaView에서 지오메트리 메시를 엽니 다. 그런 다음 유형 "상자"의 클립 (또는 슬라이스)을 사용합니다. 이렇게하면 ParaView 3D 위젯 상자를 마우스로 신속하고 대화식으로 이동, 회전 및 크기 조정할 수 있습니다.

속성 패널은 항상 나에게 보여 상자 매개 변수 :

  • 위치
  • 회전
  • 규모

지금까지 너무 좋아,하지만 난 방법을 사용할 수 있습니다 이러한 ParaView 밖에서이 상자로 작업 할 값 (예 : MATLAB)? 상자의 실제 좌표는 무엇입니까?

Ventricle Placement Below Atria

답변

0

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