Извлечение подмножества набора данных объема
[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 в пределах указывают, что объем не должен быть обрезан вдоль той оси.)
Массивы XY, и Z задайте координаты для объема V. Подобъем возвращен в NV и координаты подобъема даны в NX, NY, и NZ.
[Nx,Ny,Nz,Nv] = subvolume(V,limits) принимает массивы XY, и Z заданы как
[X,Y,Z] = meshgrid(1:N,1:M,1:P)
где [M,N,P] = size(V).
Nv = subvolume(...) возвращает только подобъем.
Загрузите mri набор данных, который является набором срезов 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