Извлеките наклонную часть из 3D объемных данных
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) в объемных данных, составляет верхний левый пиксель в плоскости изображения. Чтобы заполнить пиксельные значения в плоскости изображения, начните с угла среза и считайте значения интенсивности в порядке сканирования слева направо, сверху вниз.
