Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (Bashardoust_SystCoord) Browse

R8218/ShoulderCase/@DicomVolume/DicomVolume.m
a90c6220da75Bashardoust_SystCoord

Last Change
Raw File

DicomVolume.m
View Options

classdef (Abstract) DicomVolume < handle
properties
volume = []
spatial = []
dicomInfo = []
dicomFolderPath = []
end
methods
function set.volume(obj,volume)
assert(length(size(volume))>=3,'Volume must be a 3D array')
obj.volume = squeeze(volume);
end
function set.spatial(obj,spatial)
requiredStructureField = {'PatientPositions','PixelSpacings','PatientOrientations'};
assert(all(contains(fields(spatial),requiredStructureField)),'Invalid spatial structure');
obj.spatial = spatial;
end
function loadDataFromFolder(obj,dicomFolderPath)
assert(obj.pathContainsDicomFiles(dicomFolderPath));
obj.dicomFolderPath = dicomFolderPath;
[obj.volume,obj.spatial,~] = dicomreadVolume(dicomFolderPath);
obj.loadDicomInfoFromFolder(dicomFolderPath);
end
function loadDicomInfoFromFolder(obj,dicomFolderPath)
assert(obj.pathContainsDicomFiles(dicomFolderPath));
dicomFiles = dir(fullfile(dicomFolderPath,'*.dcm'));
obj.dicomInfo = dicominfo(fullfile(dicomFiles(1).folder,dicomFiles(1).name));
end
function output = pathContainsDicomFiles(obj,path)
output = false;
assert(isfolder(path),'Provided argument is not a valid folder path');
dicomFiles = dir(fullfile(path,'*.dcm'));
assert(not(isempty(dicomFiles)),'No dicom file found there %s',path);
output = true;
end
function output = getPointIndexInVolume(obj,point)
% the output are the three indices of the given point in the given volume.
% i,j,k: output indices
k = find(abs(obj.spatial.PatientPositions(:,3)-point(3)) == min(abs(obj.spatial.PatientPositions(:,3)-point(3))));
i = round( (point(1) - obj.spatial.PatientPositions(k,1)) / obj.spatial.PixelSpacings(k,1) );
j = round( (point(2) - obj.spatial.PatientPositions(k,2)) / obj.spatial.PixelSpacings(k,2) );
assert(all([i j k] < size(obj.volume)),'Indices found are outside the volume boundaries');
output = [i,j,k];
end
end
end
c4science · Help