Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (Bashardoust_SystCoord) Browse

R8218/ShoulderCase/@DicomVolumeNormaliser/DicomVolumeNormaliser.m
b55b6a89fc1aBashardoust_SystCoord

Last Change
Raw File

DicomVolumeNormaliser.m
View Options

classdef DicomVolumeNormaliser < handle & DicomVolume
properties
resolution = [];
normalisedVolume = [];
end
methods
function obj = DicomVolumeNormaliser(varargin)
if (nargin == 1)
obj.loadDataFromFolder(varargin{1});
end
end
function setMinimalResolution(obj)
zSpacing = obj.spatial.PatientPositions(2:end,3)-obj.spatial.PatientPositions(1:end-1,3);
minZ = min(zSpacing,[],'all');
minXY = min(obj.spatial.PixelSpacings,[],'all');
obj.resolution = min([minXY minZ],[],'all')/2;
end
function normalise(obj)
normalisedSizeZ = round(abs(obj.spatial.PatientPositions(1,3) - obj.spatial.PatientPositions(end,3))/obj.resolution);
normalisedSizeX = round(max(obj.spatial.PixelSpacings(:,1),[],'all')*size(obj.volume,1)/obj.resolution);
normalisedSizeY = round(max(obj.spatial.PixelSpacings(:,2),[],'all')*size(obj.volume,2)/obj.resolution);
obj.normalisedVolume = imresize3(obj.volume,[normalisedSizeX normalisedSizeY normalisedSizeZ]);
end
function OLDnormalise(obj)
normalisedSizeZ = round(abs(obj.spatial.PatientPositions(1,3) - obj.spatial.PatientPositions(end,3))/obj.resolution);
normalisedSizeX = round(max(obj.spatial.PixelSpacings(:,1),[],'all')*size(obj.volume,1)/obj.resolution);
normalisedSizeY = round(max(obj.spatial.PixelSpacings(:,2),[],'all')*size(obj.volume,2)/obj.resolution);
obj.normalisedVolume = zeros(normalisedSizeX,normalisedSizeY,normalisedSizeZ);
for x = 1:size(obj.normalisedVolume,1)
disp(sprintf('%d', round(100*x/size(obj.normalisedVolume,1))))
for y = 1:size(obj.normalisedVolume,2)
for z = 1:size(obj.normalisedVolume,3)
pointPosition = obj.spatial.PatientPositions(1,:)+obj.resolution*[x y z];
try
pointIndices = obj.getPointIndexInDicomVolume(pointPosition);
catch
[x y z]
pointPosition
end
obj.normalisedVolume(x,y,z) = obj.volume(pointIndices(1),pointIndices(2),pointIndices(3));
end
end
end
end
end
end
c4science · Help