Расчет геометрии торцевой крышки изо-поверхности
fvc = isocaps(X,Y,Z,V,isovalue)
fvc = isocaps(V,isovalue)
fvc = isocaps(...,'enclose')
fvc = isocaps(...,'whichplane')
[f,v,c] = isocaps(...)
isocaps(...)
fvc = isocaps(X,Y,Z,V,isovalue) вычисляет геометрию торца изо-поверхности для данных объема V при значении изо-поверхности isovalue. Множества X, Y, и Z определить координаты для объема, V.
Структура fvc содержит данные грани, вершины и цвета для конечных крышек и может быть передан непосредственно в patch команда.
fvc = isocaps(V,isovalue) предполагает массивы X, Y, и Z определяются как [X,Y,Z] = meshgrid(1:n,1:m,1:p) где [m,n,p] = size(V).
fvc = isocaps(...,' указывает, включают ли конечные ограничения значения данных выше или ниже значения, указанного в enclose') isovalue. 'enclose' параметр может быть либо 'above' (по умолчанию) или 'below'.
fvc = isocaps(...,' указывает, на каких плоскостях следует рисовать концевые крышки. Возможные значения для whichplane')'whichplane' являются 'all' (по умолчанию), 'xmin', 'xmax', 'ymin', 'ymax', 'zmin', или 'zmax'.
[f,v,c] = isocaps(...) возвращает данные грани, вершины и цвета для конечных крышек в трех массивах вместо структуры fvc.
isocaps(...) без выходных аргументов рисует фрагмент с вычисленными гранями, вершинами и цветами.
В этом примере используется набор данных, который представляет собой совокупность фрагментов МРТ человеческого черепа. Он иллюстрирует использование изокапов для прорисовки концевых колпачков на этом вырезанном объеме.
Красная изоповерхность показывает контур объема (черепа), а концевые колпачки показывают, что находится внутри объема.
Исправление, созданное на основе конечных данных (p2) использует интерполированную окраску лица, что означает, что серый цвет и источники света определяют, как он окрашен. Накладка изоповерхности (p1) использовал плоский красный цвет лица, на который влияют огни, но не использует карту цвета.
load mri D = squeeze(D); D(:,1:60,:) = []; p1 = patch(isosurface(D, 5),'FaceColor','red',... 'EdgeColor','none'); p2 = patch(isocaps(D, 5),'FaceColor','interp',... 'EdgeColor','none'); view(3) axis tight daspect([1,1,.4]) colormap(gray(100)) camlight left camlight lighting gouraud isonormals(D,p1)
isonormals | isosurface | reducepatch | reducevolume | smooth3 | subvolume