function output = getContour(obj)
  muscleMask = obj.loadMask();
  slicePixelCoordinates = obj.loadPixelCoordinates();

  output.mask = getContourMask(muscleMask);
  output.indices = getContourIndices(muscleMask);
  output.coordinates = getContourCoordinates(muscleMask,slicePixelCoordinates);
end

function output = getContourMask(muscleMask)
  contourMask = bwmorph(muscleMask,"remove");
  output = contourMask;
end

function output = getContourCoordinates(muscleMask,pixelCoordinates)
  contourIndices = getContourIndices(muscleMask);
  contourCoordinates = [];
  for row = 1:size(contourIndices,1)
    contourCoordinates = [contourCoordinates; pixelCoordinates{contourIndices(row,1),contourIndices(row,2)}];
  end
  output = contourCoordinates;
end

function output = getContourIndices(muscleMask)
  contourMask = getContourMask(muscleMask);
  [row, column] = ind2sub(size(contourMask),find(contourMask));
  contourIndices = [column row];
  output = contourIndices;
end