Извлеките наклонный срез из 3-D объемных данных
B = obliqueslice(V,point,normal)V. Срез извлечен со ссылкой на данную точку на объеме и векторе нормали. Режущая плоскость перпендикулярна вектору нормали и проходит через заданную точку.
Для получения информации о том, как срез извлечен относительно данной точки и нормального, смотрите Наклонное Разрезание. Ориентация извлеченного среза в плоскости изображения зависит от своего положения в 3-D координатном пространстве. Для получения дополнительной информации смотрите Ориентацию Среза в Плоскости изображения.
B = obliqueslice(___,Name,Value)
Загрузите 3-D объемный набор данных в рабочую область.
load mriУдалите одноэлементные размерности при помощи squeeze функция.
V = squeeze(D);
Отобразите горизонтальные срезы данных при помощи montage функция.
montage(V,map,'Size',[3 9]);
Задайте точку в объеме для среза, чтобы пройти.
point = [73 50 15.5];
Задайте вектор нормали в 3-D координатном пространстве.
normal = [0 15 20];
Извлеките срез из объемных данных. Срез перпендикулярен вектору нормали и проходит через заданную точку.
[B,x,y,z] = obliqueslice(V,point,normal);
Отобразите извлеченный срез в 3-D координатном пространстве.
figure surf(x,y,z,B,'EdgeColor','None','HandleVisibility','off'); grid on view([-38 12]) colormap(gray) xlabel('x-axis') ylabel('y-axis'); zlabel('z-axis'); title('Slice in 3-D Coordinate Space')
Постройте точку и вектор нормали.
hold on plot3(point(1),point(2),point(3),'or','MarkerFaceColor','r'); plot3(point(1)+[0 normal(1)],point(2)+[0 normal(2)],point(3)+[0 normal(3)], ... '-b','MarkerFaceColor','b'); hold off legend('Point in the volume','Normal vector')
Отобразите извлеченный срез в плоскости изображения.
figure
imshow(B,[])
title('Slice in Image Plane')
Загрузите 3-D объемный набор данных в рабочую область.
s = load(fullfile(toolboxdir('images'),'imdata','BrainMRILabeled','images','vol_001.mat')); V = s.vol;
Отобразите горизонтальные срезы данных при помощи montage функция.
montage(V,'Indices',12:118,'Size',[8 12],'DisplayRange',[]);
Задайте вектор нормали к плоскости в 3-D координатном пространстве.
normal = [20 0 10];
Извлеките несколько срезов вдоль направления вектора нормали с помощью for цикл. В каждой итерации:
Задайте точку, что срез должен пройти.
Извлеките срез, задав выходной размер к 'Full' и значение заливки для дополнения пикселей как 255. Извлеченные срезы перпендикулярны вектору нормали и проходят через заданную точку.
Отобразите извлеченные срезы.
sliceIdx = 10:5:180; figure for s = 1:length(sliceIdx) pt = [sliceIdx(s) 150 80]; [B,x,y,z] = obliqueslice(V,pt,normal,'OutputSize','Full','FillValues',255); Bslices(:,:,s) = B; % Display the slice in 3-D coordinate space subplot('Position',[0.11 0.36 0.38 0.5]) surf(x,y,z,B,'EdgeColor','None','HandleVisibility','off'); grid on view([-24 12]) colormap(gray) xlabel('x-axis') ylabel('y-axis'); zlabel('z-axis'); zlim([0 155]); ylim([0 250]); xlim([0 250]); title('Slice in 3-D Coordinate Space') % Plot the point and the normal vector. hold on plot3(pt(1),pt(2),pt(3),'or','MarkerFaceColor','r') plot3( ... pt(1)+[-normal(1) normal(1)], ... pt(2)+[-normal(2) normal(2)], ... pt(3)+[-normal(3) normal(3)], ... '-b','MarkerFaceColor','b') legend('Point in the volume','Normal vector','Position',[0.1 0.12 0.3 0.08]) hold off % Display the extracted slice. subplot('Position',[0.6 0.37 0.34 0.49]) imshow(B,[]) title('Slice in Image Plane') pause(0.5); end
Отобразите извлеченные фрагменты изображения при помощи montage функция.
figure montage(Bslices,'Size',[5 7],'DisplayRange',[]);
Учитывая точку (px, py, pz) и вектор нормали (a, b, c), функция решает плоское уравнение a (x-px) +b (y-py) +c (z-pz) = 0
Точка (px, py, pz) находится в объемных данных. Режущая плоскость перпендикулярна вектору нормали и проходит через данную точку.
Ориентация извлеченного среза в плоскости изображения зависит от своего наклонного угла относительно горизонтальных и вертикальных плоскостей.
Угол среза, который находится близкий началу координат из последнего среза, (0, 0, P), в объемных данных составляет верхний левый пиксель в плоскости изображения. Чтобы заполнить пиксельные значения в плоскости изображения, запустите с угла среза и считайте значения интенсивности в слева направо, порядок скана от начала до конца.
