Вычисление геометрии заглушки изоповерхности
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.
Struct 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(...) возвращает поверхность, вершину и цветные данные для заглушек в трех массивах вместо struct fvc.
isocaps(...) без выходных аргументов чертит закрашенную фигуру с вычисленными поверхностями, вершинами и цветами.
Этот пример использует набор данных, который является набором срезов MRI человеческого черепа. Это иллюстрирует использование изопрописных букв, чтобы чертить заглушки на этом срезанном объеме.
Красная изоповерхность показывает схему объема (череп), и заглушки показывают то, что в объеме.
Закрашенная фигура, созданная из данных заглушки (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