Model reconstruction from CT slices
10 ビュー (過去 30 日間)
古いコメントを表示
Hello everyone,
I apologize if this topic was repeated but I have not found the solution yet. I have 537 DICOM images, and i can extract the coordinates (x,y) of the bony region (the intensity which is equal to 1). Now i need reconstruct all images in a 3D model and i have no idea how to do it. I need pair/match the pixels that are 1?
Best regards and thanks for your attention,
Nuno
0 件のコメント
採用された回答
Christiaan
2015 年 3 月 11 日
編集済み: Christiaan
2015 年 3 月 11 日
Dear Nuno,
You can load and plot the individual DICOM images with the hand made code below. Make sure that the path of the m-file is equal to the path of the images. As a test I used the images from this thread .
clc;clear all;close all;
% Read the series of images.
for p=1:9
filename = sprintf('00000%d.dcm', p);
X(:,:,1,p) = dicomread(filename);
X(:,:,1,p) = imadjust(X(:,:,1,p),[0 0.022],[]); % increase contrast
end
for p=10:42
filename = sprintf('0000%d.dcm', p);
X(:,:,1,p) = dicomread(filename);
X(:,:,1,p) = imadjust(X(:,:,1,p),[0 0.022],[]); % increase contrast
end
numberOfSlices = p;
figure(1); montage(X,[]);
% to plot for example one image (here nr. 10)
figure(2); imshow(X(:,:,1,10))
The second step would be to stack the 2D images together. In this thread multiple possibilities are gives how to do this. A nice feature is sliding viewer. If you download that m-file and put the files in the map where also the images are located and then write in your matlab-prompt:
D=double(squeeze(X));
slidingviewer(D);
you will see a 3D image of your CT scan.
Good luck! Christiaan
その他の回答 (1 件)
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!