function [sphereValue image recoveryCoeff] = getSpheres(info, image)
    imageSize= size(image);
    %pixelSpacing = info.PixelSpacing(1);
    pixelSpacing = 3.1819;
    %matrixSize = double(info.Columns(1)/2);
    
    %create the radius and center of the spheres
    %this is dependent of the image we have
    %it changes with pixel- and matrix size
    %later on it has to be automatised
    radius37 = 18.5/3.1819;
    %center37 = [matrixSize+57+2/pixelSpacing,matrixSize+2, 44];
    %center37=[(matrixSize-28)/pixelSpacing, matrixSize-15.4 ,44];
    center37 = [130.5,149,37];
    
    radius13 = 6.5/pixelSpacing;
    %center17 = [matrixSize-57/pixelSpacing, matrixSize, 44];
    center13=[145,121,37];

    radius28 = 14/pixelSpacing;
    %center28 = [ matrixSize+28/pixelSpacing, matrixSize+15.4 ,44];
    center28=[114.5, 140, 37];
    
    radius17 = 8.5/pixelSpacing;
    %center22 = [matrixSize-28/pixelSpacing, matrixSize+15.4 ,44];
    center17=[129,112,37];
    
    radius22 = 11/pixelSpacing;
    %center13 = [matrixSize+57+2/pixelSpacing,matrixSize+2, 44];
    center22=[113.5, 123, 37];

    radius10 = 5/pixelSpacing;
    %center10 = [matrixSize+28+2/pixelSpacing, matrixSize-15.4+2 ,44];
    center10=[145,139,37];
    % Create a logical image of a circle with specified
    % diameter, center, and image size.
    % First create the image.
    imageSizeX = imageSize(1);
    imageSizeY = imageSize(2);
    imageSizeZ = imageSize(3);
    [columnsInImage rowsInImage depth] = meshgrid(1:imageSizeX, 1:imageSizeY, 1:imageSizeZ);
    %in z-direction I have a different pixel spacing as in x- and
    %y-direction
    %for this I need this correction factor
    correctionZ=(2/3.1819)*(2/3.1819);
    %create the spheres
    circlePixels37 = (rowsInImage - center37(2)).^2 ...
                 + (columnsInImage - center37(1)).^2 ...
                 + correctionZ*(depth - center37(3)).^2 <= radius37.^2;
    circlePixels17 = (rowsInImage - center17(2)).^2 ...
               + (columnsInImage - center17(1)).^2 ...
                 + correctionZ*(depth - center17(3)).^2 <= radius17.^2;
    circlePixels22 = (rowsInImage - center22(2)).^2 ...
                 + (columnsInImage - center22(1)).^2 ...
                 + correctionZ*(depth - center22(3)).^2 <= radius22.^2;

    circlePixels28 = (rowsInImage - center28(2)).^2 ...
                 + (columnsInImage - center28(1)).^2 ...
                 + correctionZ*(depth - center28(3)).^2 <= radius28.^2;        

    circlePixels13 = (rowsInImage - center13(2)).^2 ...
               + (columnsInImage - center13(1)).^2 ...
               + correctionZ*(depth - center13(3)).^2 <= radius13.^2; 

    circlePixels10 = (rowsInImage - center10(2)).^2 ...
               + (columnsInImage - center10(1)).^2 ...
               + correctionZ*(depth - center10(3)).^2 <= radius10.^2; 
    %get the values inside the spheres     
    sphereValue=image(circlePixels37| circlePixels17 | circlePixels22 | circlePixels28 | circlePixels13 | circlePixels10);
      
     sphereValue37=image(circlePixels37);
     sphereValue17=image(circlePixels17);  
     sphereValue22=image(circlePixels22);  
     sphereValue28=image(circlePixels28);
     sphereValue13=image(circlePixels13)  
     sphereValue10=image(circlePixels10);  
     recoveryCoeff37 = double(max(double(sphereValue37))/double(23000.0));
     recoveryCoeff17 = double(max(double(sphereValue17))/double(23000.0));
     recoveryCoeff22 = double(max(double(sphereValue22))/double(23000.0));
     recoveryCoeff28 = double(max(double(sphereValue28))/double(23000.0));
     recoveryCoeff13 = double(max(double(sphereValue13))/double(23000.0));
     recoveryCoeff10 = double(max(double(sphereValue10))/double(23000.0));
%      recoveryCoeff37 = mean(sphereValue37(sphereValue37>5000))/23000;
%      recoveryCoeff17 = mean(sphereValue17(sphereValue17>5000))/23000;
%      recoveryCoeff22 = mean(sphereValue22(sphereValue22>5000))/23000;
%      recoveryCoeff28 = mean(sphereValue28(sphereValue28>5000))/23000;
%      recoveryCoeff13 = mean(sphereValue13(sphereValue13>5000))/23000;
%      recoveryCoeff10 = mean(sphereValue10(sphereValue10>5000))/23000;
    %set the spheres to 0 in the image so that a calculation of the
    %background values becomes more easy
        
    image(circlePixels37| circlePixels17 | circlePixels22 | circlePixels28 | circlePixels13 | circlePixels10) = 0;
%     figure()    
%     h =imshow(image(:,:,37)*2)
     %g =imshow(image(:,:,34)*2)
%     savefig(h, 'D:/EvaluationSimulation_buggy/Simulation_1_10/TRUEXM256/120s/example.png');

    recoveryCoeff= zeros(1,6, 'double');
    recoveryCoeff(1,1) = recoveryCoeff37;
    recoveryCoeff(1,2) = recoveryCoeff28;
    recoveryCoeff(1,3) = recoveryCoeff22;
    recoveryCoeff(1,4) = recoveryCoeff17;
    recoveryCoeff(1,5) = recoveryCoeff13;
    recoveryCoeff(1,6) = recoveryCoeff10;
    recoveryCoeff
end