Filter löschen
Filter löschen

Info

Diese Frage ist geschlossen. Öffnen Sie sie erneut, um sie zu bearbeiten oder zu beantworten.

how to form volume from the 2d slice of gated spect images

2 Ansichten (letzte 30 Tage)
andhavarapu lokesh
andhavarapu lokesh am 24 Jan. 2017
Geschlossen: MATLAB Answer Bot am 20 Aug. 2021
dear members iam new to matlab want to form a volume from 2d slices where each dcm file has 72 frames where 1-35 one volume and 36 to 72 other volume is there any way to form volume reconstruction or rendering
if true
projectdir = 'E:\SHIVA BACKUP\THYROID\P1\newcodes\data1\13002';
X = zeros(128, 128, 1, 72,'uint16');
% Read the series of images.
% for p=1
p=1;
thisfile = sprintf('IM_%d.dcm', p);
filename = fullfile( projectdir, thisfile );
imdata = dicomread(filename);
imsize = size(imdata);
if ~isequal( imsize, [128 128 1 72] )
fprintf('file is unexpected size %s instead of [128 128 1 72], skipping "%s"\n', mat2str(imsize), filename);
else
X(:, :, :, :,p) = imdata;
end
figure();
montage(X(:,:,:,:,p), [0 100])
s=[];
v=zeros(128,128,72);
sum=0;
v=squeeze(X);
for i=1:1:35
v(:,:,i)=X(:,:,:,i);
figure()
imshow(v(:,:,i),[0 100])
end
volume=sum*8;
isoval=15;
fig = figure();
ax = axes('Parent', fig);
fv=isosurface(v,isoval);
hiso = patch(fv, ...
'Parent', ax, 'FaceAlpha', 0.74, ...
'FaceColor', [0.75,0.75,0.75],'EdgeColor','red');
lighting(ax, 'phong');
%lightangle(ax, 90, 30);
rotate3d(ax, 'on');
title( sprintf('p = %d', p) )
[K,V] = convhull(fv.vertices);
end

Antworten (0)

Diese Frage ist geschlossen.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by