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