ParaView Clip/Slice, Box Параметры: Значение положения, Поворот, Масштаб; Реконструировать коробку в Matlab

Question

Я хочу определить 3D-бокс определенного размера и местоположения рядом с заданной триангулированной геометрией.

Поэтому я открываю геометрическую сетку в ParaView. Затем я использую клип (или фрагмент) типа "Box". Таким образом, я получаю виджет ParaView 3D для блока, который можно перемещать, поворачивать и изменять размер быстро и в интерактивном режиме с помощью мыши.

Панель свойств всегда показывает мне параметры окна:

• Позиция
• вращение
• Масштаб

Пока все хорошо, но как я могу использовать эти значения для работы с этим полем вне ParaView, например в MATLAB? Каковы фактические координаты коробки?

Answer 1

Мы копируем свойства блока из 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.

% 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)