Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (master) Browse

R8218/ShoulderCase/@MaskSubdividerWithSquareGrid/MaskSubdividerWithSquareGrid.m
a850000c9cd2master

Last Change
Raw File

MaskSubdividerWithSquareGrid.m
View Options

classdef MaskSubdividerWithSquareGrid < handle
properties
mask = [];
maskLimit = [];
xGrid = [];
xResolution = 1;
yGrid = [];
yResolution = 1;
end
methods
function obj = MaskSubdividerWithSquareGrid(mask)
obj.mask = mask;
[obj.xGrid, obj.yGrid] = meshgrid(1:size(mask, 1), 1:size(mask, 2));
[maskPixelY, maskPixelX] = find(mask);
obj.maskLimit.top = min(maskPixelY);
obj.maskLimit.bottom = max(maskPixelY);
obj.maskLimit.left = min(maskPixelX);
obj.maskLimit.right = max(maskPixelX);
end
function setXResolutionInPixels(obj, resolutionInPixels)
resolutionInPixels = round(resolutionInPixels);
obj.xResolution =...
min(obj.maskLimit.right - obj.maskLimit.left,...
max(1, resolutionInPixels));
end
function setXResolutionInMm(obj, resolutionInMm, pixelSizeX)
obj.setXResolutionInPixels(resolutionInMm/pixelSizeX);
end
function setYResolutionInPixels(obj, resolutionInPixels)
resolutionInPixels = round(resolutionInPixels);
obj.yResolution =...
min(obj.maskLimit.bottom - obj.maskLimit.top,...
max(1, resolutionInPixels));
end
function setYResolutionInMm(obj, resolutionInMm, pixelSizeY)
obj.setYResolutionInPixels(resolutionInMm/pixelSizeY);
end
function output = getDivision(obj, top, bottom, left, right)
division =...
(obj.xGrid >= left)...
& (obj.xGrid <= right)...
& (obj.yGrid >= top)...
& (obj.yGrid <= bottom);
output = division;
end
function output = getMaskSubdivisions(obj)
maskSubdivisions = [];
for x = obj.maskLimit.left:obj.xResolution:obj.maskLimit.right
for y = obj.maskLimit.top:obj.yResolution:obj.maskLimit.bottom
maskSubdivisions(:,:,end+1) =...
obj.mask...
& obj.getDivision(...
y, (y + obj.yResolution-1), x, (x + obj.xResolution-1));
end
end
output = obj.removeEmptySubdivisions(maskSubdivisions);
end
function output = removeEmptySubdivisions(obj, subdivisions)
subdivisionsToRemove = [];
for i = 1:size(subdivisions, 3)
if not(any(subdivisions(:,:,i), "all"))
subdivisionsToRemove(end+1) = i;
end
end
subdivisions(:,:,subdivisionsToRemove) = [];
output = subdivisions;
end
end
end
c4science · Help