Уменьшайте число элементов в наборе данных объема
[nx,ny,nz,nv] = reducevolume(X,Y,Z,V,[Rx,Ry,Rz])
[nx,ny,nz,nv] = reducevolume(V,[Rx,Ry,Rz])
nv = reducevolume(...)
[nx,ny,nz,nv] = reducevolume(X,Y,Z,V,[Rx,Ry,Rz]) уменьшает число элементов в объеме путем сохранения каждого Rxth элемент в направлении X, каждом Ryth элемент в направлении Y и каждый Rzth элемент в z направлении. Если скалярный R используется, чтобы указать на сумму или сокращение вместо трехэлементного вектора, MATLAB® программное обеспечение принимает сокращение, чтобы быть [R R R].
Массивы XY, и Z задайте координаты для объема V. Уменьшаемый объем возвращен в nv, и координаты уменьшаемого объема возвращены в nx, ny, и nz.
[nx,ny,nz,nv] = reducevolume(V,[Rx,Ry,Rz]) принимает массивы XY, и Z заданы как [X,Y,Z] = meshgrid(1:n,1:m,1:p), где [m,n,p] = size(V).
nv = reducevolume(...) возвращает только уменьшаемый объем.
Этот пример использует набор данных, который является набором срезов MRI человеческого черепа. Эти данные обрабатываются во множестве путей:
4-D массив сжат (squeeze) в три измерения и затем уменьшаемый (reducevolume) так, чтобы то, что остается, было каждым четвертым элементом в направлениях X и Y и каждым элементом в z направлении.
Уменьшаемые данные сглаживаются (smooth3).
Схема черепа является изоповерхностью, сгенерированной как закрашенная фигура (p1) чьи нормали вершин повторно вычисляются, чтобы улучшить внешний вид, когда подсветка применяется (patchизоповерхностьизонормали).
Вторая закрашенная фигура (p2) с интерполированной поверхностью цвет чертит заглушки (FaceColorизопрописные буквы.
Представление объекта установлено (viewосьdaspect).
Полутоновая палитра с 100 элементами обеспечивает окраску для заглушек (colormap).
Добавление света справа от камеры освещает объект (camlightподсветка).
load mri D = squeeze(D); [x,y,z,D] = reducevolume(D,[4,4,1]); D = smooth3(D); 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,.4]) colormap(gray(100)) camlight lighting gouraud
isosurface | isocaps | isonormals | smooth3 | subvolume | reducepatch