Извлечение подмножества набора данных объема
[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] (Любой NaN s в пределах указывают, что объем не должен быть обрезан вдоль той оси.)
Массивы 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, который является набором срезов 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