Извлечь подмножество набора данных тома
[Nx,Ny,Nz,Nv] = subvolume(X,Y,Z,V,limits)
[Nx,Ny,Nz,Nv] = subvolume(V,limits)
Nv = subvolume(...)
[Nx,Ny,Nz,Nv] = subvolume(X,Y,Z,V,limits) извлекает подмножество набора данных тома V с использованием указанной выровненной оси limits. limits = [xmin,xmax,ymin, ymax,zmin,zmax] Любой NaNs в пределах указывают, что объем не должен обрезаться вдоль этой оси.)
Множества X, Y, и Z определить координаты для объема, V. Объем возвращается в NV и координаты подтома приведены в NX, NY, и NZ.
[Nx,Ny,Nz,Nv] = subvolume(V,limits) предполагает массивы X, Y, и Z определяются как
[X,Y,Z] = meshgrid(1:N,1:M,1:P)
где [M,N,P] = size(V).
Nv = subvolume(...) возвращает только объем.
Загрузить mri набор данных, который представляет собой коллекцию срезов МРТ человеческого черепа. Затем:
Выдавливание четырехмерного массива в три измерения и извлечение подмножества с помощью subvolume.
Используйте пластырь для отображения изоповерхности контура черепа.
Добавьте второй фрагмент с интерполированными цветами граней.
Задайте вид объекта.
Измените карту цветов на карту цветов в оттенках серого.
Добавьте источники света справа и слева от камеры для освещения объекта.
load mri D = squeeze(D); [x,y,z,D] = subvolume(D,[60,80,nan,80,nan,nan]); p1 = patch(isosurface(x,y,z,D, 5),... 'FaceColor','red','EdgeColor','none'); isonormals(x,y,z,D,p1); p2 = patch(isocaps(x,y,z,D, 5),... 'FaceColor','interp','EdgeColor','none'); view(3); axis tight; daspect([1 1 0.4]) colormap(gray(100)) camlight right; camlight left; lighting gouraud
isocaps | isonormals | isosurface | reducepatch | reducevolume | smooth3