function output = getPointsIndicesY(obj, pointsCoordinatesY)
  indicesY = round(...
    (pointsCoordinatesY - obj.spatial.PatientPositions(1,2))...
    / obj.spatial.PixelSpacings(1,2));
  output = indicesY;
end