Segmentation of Collarbone in CT - volume

3 次查看(过去 30 天)
Alba
Alba 2020-5-29
评论: Alba 2020-6-2
Hello,
I'm new to matlab and my task is to segment the trachea (wind pipe) and the clavicle bones (Collarbones) in different CT volumes. I could finish the segmentation of trachea with success. In the second part, I could extract the entire set of bone within the volume and I have stuck at there. I couldn't extract both collarbones from the rest of bones. As I'm supposed to measue the lenght of at lease one of these bones, could you please enlight me on how to segment these two bones?
The bone extaction code I wote so far is,
load 3D_Configs.mat; % Load the Configs
% Get the folder path
folder = ('./Set_No_2');
d = dir(folder + "/*.dcm");
slice_location = [];
for k = 1:length(d)
info = dicominfo(folder + "/" + d(k).name);
slice_location(k) = info.SliceLocation;
end
% Read the DICOM volume
[V,sp] = dicomreadVolume(folder);
V = squeeze(V);
slice_location = sort(slice_location);
scale_factors = [sp.PixelSpacings(1,:) slice_location(2)-slice_location(1)]; % Rescale the volume based on the DICOM metadata
%volumeViewer(V,'ScaleFactors',scale_factors)
% Visualise the skeleton volume
vs = volshow(V,config,'ScaleFactors',scale_factors);
(I tried creating various volume masks with the help of Image Segmenter and tried to segment the collarbones. But all failed.)
(The "3D_Configs.mat" file is a config extract saved from the volumeViewer with cam angle. Attached here.).
The "Set_No_2" is the DICOM image folder.
Code developed to segment the trachea in 2D slices of DICOM images is also attached here.
Thanks a lot in advance.
Alba

回答(1 个)

Dinesh Yadav
Dinesh Yadav 2020-6-2
Hi Alba,
I am aware of one method i.e when you know the location of points for the region of interest (i.e. collar bone) you can extract those points using imcrop3.
There are two more cropping methods randomCropWindow3d and centreCropWindow3d but using this will take trail and error if we dont know the co-ordinates of the region of interest.
Hope it helps.

类别

Help CenterFile Exchange 中查找有关 3-D Volumetric Image Processing 的更多信息

产品


版本

R2020a

Community Treasure Hunt

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

Start Hunting!

Translated by