function calculateIntensityFeatures(pixelCoords, images, boxSize)

range = (boxSize-1)/2;

fid = fopen ('intensity_features.txt', 'a');

fprintf(fid, '%s,%s,%s,%s,%s \n','intense', 'maxInt', 'minInt', 'meanInt', 'SDInt');

%intensityFeatures = zeros( size(pixelCoords, 2) , 5);

for i=1:size(images,2)

    if i==1 || (strcmp(images(i-1), images(i))==0)

        %CTimage = retrieveImage(char(images(i)));

        CTimage = preProcessDicom(char(images(i)));

    end

    bBox = double(CTimage(pixelCoords(2,i)-range:pixelCoords(2,i)+range,pixelCoords(1,i)-range:pixelCoords(1,i)+range));

    intense = bBox(range+1, range+1);

    maxInt = max(max(bBox));

    minInt = min(min(bBox));

    meanInt = mean2(bBox);

    SDInt = std2(bBox);

    fprintf(fid, '%d,%d,%d,%d,%d \n', intense, maxInt, minInt, meanInt, SDInt);

    %intensityFeatures(i,:) = [intense maxInt minInt meanInt SDInt];

    disp(i)

end

fclose(fid);

%xlswrite(filename, [{'intense'} {'maxInt'} {'minInt'} {'meanInt'} {'SDInt'}], 'A1:E1');

%success = xlswrite(filename, intensityFeatures, 'Sheet1', 'A2');