% spine of the mathematical phantom
function image = setSpineToZero( image )
    imageSize=size(image);
    centerSpine1_10 = [130,130];
    %centerSpine1_5 = [128,130];
    %centerSpine = [128,128];
    radiusSpine = 8;
    imageSizeX = imageSize(1);
    imageSizeY = imageSize(2);
    
    [columnsInImage rowsInImage] = meshgrid(1:imageSizeX, 1:imageSizeY);
    circlePixelsSpine = (rowsInImage - centerSpine1_10(2)).^2 ...
                 + (columnsInImage - centerSpine1_10(1)).^2  <= radiusSpine.^2;
    image(circlePixelsSpine) = 0;


end

